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PREFACE 


This  work  was  performed  for  the  Eustis  Directorate,  U.  S.  Army  Air  Mobility 
Research  and  Development  Laboratory,  Fort  Eustis,  Virginia,  under  Contract 
DAAJ02-72-C-0105 ,  Project  1F162208AA82 .  Technical  monitor  for  the  contract 
was  Mr.  Paul  Mirick.  The  testing  of  the  flexured  airfoil  model  and  re¬ 
duction  of  the  experimental  data  were  performed  at  the  United  Aircraft 
Research  Laboratories,  East  Hartford,  Connecticut.  The  airfoil  models  and 
experimental  hardware  were  designed  and  constructed  by  Mr.  W.  Hayden.  The 
testing  was  conducted  by  Mr.  P.  G.  Vogt.  Processing  of  the  stall  flutter 
data  and  correlation  with  the  UAEL  unsteady  aerodynamic  analysis  were 
performed  by  Mr.  G.  L.  Commerford.  The  investigation  of  the  effects  of 
blade  parameters  on  helicopter  stall  boundaries  was  conducted  at  the 
Sikorsky  Aircraft  Division  of  United  Aircraft  Corporation,  Stratford, 
Connecticut,  by  Mr.  R.  H.  Blackwell. 
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INTRODUCTION 


In  the  preliminary  c.esign  of  hell  copter  ■  rotors  ,  the  rotor  and  blade  proper¬ 
ties  are  established  by  applying  a  set  of  design  criteria.  Important  cri¬ 
teria  that  must  be  considered  are  those  imposed  by  retreating  blade  stall 
in  forward  flight.  The  rotor's  operating  capability  tends  to  be  determined 
by  the  degree  to  which  stall  can  be  penetrated  before  any  adverse  effects 
become  intolerable.  As  a  result,  the  form  which  the  stall-induced  limita¬ 
tion  assumes  depends  upon  the  specific  rotor  design.  For  example,  one 
rotor  may  experience  stall-ind- .ced  drag  rise  or  loss  in  lifting  capability 
or  control  power.  Another  may  be  limited  by  a  buildup  in  control  loads 
due  to  stall  flutter.  Still  others  may  experience  flatwise  or  edgewise 
stress  limits.  Most  experience  in  evolving  preliminary  design  criteria 
accounting  for  retreating  blade  stall  has  been  based  on  simole  static 
aerodynamic  criteria  such  as  maximum  angle  of  attack,  profile  drag  diver¬ 
gence,  etc.  Based  on  correlation  with  their  previous  flight  test  exper¬ 
ience,  different  organizations  have  relied  on  the  correspondence  of  their 
primary  stall  limitation  (excessive  blade  stress,  sharp  rise  in  control 
loads,  loss  of  control  power,  etc.)  with  a  particular  value  of  the  chosen 
static  aerodynamic  criterion.  Clearly,  the  reliability  of  this  type  of 
prediction  method  decreases  when  new  configurations  and/or  flight  regimes 
are  being  considered.  An  ideal  preliminary  design  method  should  predict 
the  effects  of  all  the  governing  design  parameters  on  the  various  forms 
of  stall  limitation.  Obviously,  the  fewer  the  generalized  parameters 
involved,  the  simpler  and  more  convenient  would  be  the  design  criteria. 
Progress  toward  this  goal  has  been  hampered  by  a  lack  of  understanding 
of  the  complex  rotor  aeroelastic  response  phenomena  in  stall .  Only 
recently,  with  the  incorporation  of  important  unsteady -  aerodynamic  effects, 
have  blade  response  analyses  become  capable  of  modeling  rotor  aeroelastic 
response  in  stall  (Reference  l). 

In  light  of  these  recent  advances,  it  was  appropriate  go  conduct  the  work 
described  below.  The  principal  objectives  of  this  effort  were  (l)  to  em¬ 
ploy  the  latest  available  analytical  techniques  to  define  systematically 
retreating  blade  stall  boundaries  for  a  typical  rotor,  (2)  to  assess  the 
impact  of  blade  design  parameters  on  these  boundaries,  and  (3)  where 
possible,  to  provide  a  set  of  more  generalized  stall  boundaries  useful  for 
preliminary  design  purposes.  A  further  objective  of  the  work  was  to  con¬ 
firm  the  general  accuracy  of  the  analysis  and  point  out  any  problem  areas 
through  correlation  with  flight  test  data  for  two  helicopters  and  with 
stall  flutter  data  obtained  in  a  two-dimensional  airfoil  wind  tunnel  test. 
The  correlation  of  two-dimensional  airfoil  stall  flutter  results  provided 
a  unique  opportunity  to  evaluate  different  unsteady  stall  aerodynamic 
models  and  to  gain  insight  into  stall-induced  torsional  response  phenomena 
without  the  complicating  effects  of  nonuniform  rotor  inflow  and  blade 
flexibility . 


TWO-DIMENSIONAL  DYNAMIC  STALL  EXPERIMENTS 


One  of  the  primary  objectives  of  this  investigation  was  to  provide  data 
useful  in  evaluating  the  unsteady  aerodynamic  methods  presently  used  in 
blade  aeroelastic  response  analyses.  Previous  evaluations  of  unsteady 
aerodynamic  stall  models  have  centered  on  the  correlation  between  measured 
and  calculated  lift  and  pitching  moment  hysteresis  loops  of  rigid  airfoils 
forced  through  sinusoidal  angle  of  attack  changes.  The  best  measure  of  an 
aerodynamic  method  is  its  ability  to  predict  the  stall-induced  torsional 
oscillations  experienced  on  a  full-scale  helicopter  in  forward  flight. 
However,  the  complexity  of  the  rotor  operating  environment  (blade  bending 
and  plunging  motion,  inflow  fluctuation  and  three-dimensional  effects) 
prevents  fundamental  study  of  stall  flutter  response  phenomena  on  a 
full-scale  rotor.  As  an  intermediate  approach,  the  test  program  described 
below  was  conducted.  Blade  torsional  response  to  stall  was  tested  with  a 
flexured  two-dimensional  airfoil  dynamically  scaled  to  simulate  tne  first 
torsional  mode  of  a  rotor  blade.  In  this  way  the  analytical  capability  for 
predicting  stall  flutter  response  could  be  e-  amined  without  the  uncertain¬ 
ties  of  rotor  blade  translational  motion  and  inflow.  Additionally,  this 
type  of  test  provided  the  opportunity  for  experimentally  determining  the 
effect  of  blade  section  design  parameters  on  torsional  response  behavior. 

TEST  FACILITY  AND  MODELS 


The  facility  chosen  for  the  two-dimensional  experiments  was  the  UARL 
Acoustic  Research  Wind  Tunnel.  Although  this  facility  was  specifically 
designed  for  aerodynamic  noise  research,  it  has  many  features  which  also 
made  it  useful  for  this  study.  It  is  a  low  turbulence  level,  open  circuit, 
open  jet  wind  tunnel  whose  test  section  is  surrounded  by  a  large  anechoic 
chamber.  The  test  section  of  the  tunnel  is  shown  in  figure  1.  Air  is 
drawn  from  the  atmosphere  through  a  large  bellmouth  and  then  through  the 
test  section  by  a  centrifugal  fan  at  the  exhaust  end. 

The  facility  permits  the  use  cf  models  of  adequate  size  to  obtain  Reynolds 
number  of  order  10  .  The  size  of  the  facility  permitted  airfoil  effective 
aspect  ratios  of  about  1 6  as  compared  to  an  aspect  ratio  of  9-5  in  the 
tests  of  Reference  1. 

The  airfoil  models  were  required  to  be  as  stiff  as  possible  along  their 
span  in  order  to  maintain  the  two-dimensionality  of  the  test.  A  second 
requirement  was  that  they  be  light  weight  so  that  they  would  be  reasonably 
representative  of  scaled  helicopter  blade  mass  and  inertia  properties.  To 
satisfy  these  requirements,  models  were  constructed  with  fiberglass  skins 
over  a  hollow  steel  spar  with  end  ribs,  and  foam-type  filler  material.  A 
schematic  of  the  model  construction  is  shown  in  Figure  2.  Two  different 
airfoil  sections,  one  an  NACA  0012  and  one  an  SC  1095  >  were  fabricated. 

In  each  case  the  airfoil  chord  was  0.5  ft  and  the  span  was  1.75  fn.  Max¬ 
imum  deviation  between  designed  and  measured  airfoil  thicMess  was  .005  c. 
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Each  end  of  the  airfoil  had  a  drive  shaft  which  could  be  moved  to  place 
the  pitch  axis  at  two  chordwise  locations  (.22c  and  .25c).  The  airfoil 
was  driven  from  one  side  of  the  wind  tunnel  by  a  variable  length  torque 
bar  which  served  as  the  torsional  spring  used  to  set  the  blade  pitch 
frequency.  The  opposite  side  of  the  airfoil  had  a  larger  shaft  mounted 
in  a  pair  of  roller  bearings  in  tiie  tunnel  side  wall.  It  extended  beyond 
the  outer  edge  of  the  wind  tunnel  and  had  the  rotor  of  a  capacitive  posi¬ 
tion  measuring  device  affixed  to  its  end.  The  capacitive  angular  trans¬ 
ducer  was  described  in  detail  in  Reference  1.  The  transducer  is  a 
parallel  plate  capacitor  of  concentric  cyclindric  configuration  whose  out¬ 
put  signal  is  proportional  to  the  overlap  area  (Figure  3).  The  stator  or 
fixed  portion  of  this  transducer  was  attached  to  the  outside  of  the  wind 
tunnel  side  wall.  The  free  shaft,  end  was  threaded  so  that  an  inertial 
disk  couj.d  be  mounted  on  the  shaft  to  increase  airfoil  torsional  inertia 
by  approximately  50%. 

The  oscillatory  mechanism  which  was  described  in  Reference  1  is  shown 
schematically  in  Figure  U  along  with  the  details  of  the  transmission  sys¬ 
tem.  The  spring  connector  between  the  airfoil  and  the  transmission  was 
modified  for  this  program  as  shown  in  the  cutaway  drawing  in  Figure  5. 

The  drive  shaft  from  the  transmission  was  a  thick -walled  hollow  tube.  A 
collet  mechanism  was  mounted  inside  this  Lube  in  a  kevvay  .  "he  collet 
held  the  smaller  diameter  shaft  (tne  torque  bar)  of  the  airfoil  and  coula 
be  moved  axially  to  clamp  this  shaft.  The  far  end  of  the  hollow  shaft 
was  supported  in  a  large  bearing  in  the  tunnel  side  wall. 

The  mean  angle  of  attack  of  the  airfoil  was  set  by  adjusting  the  relative 
positions  of  the  solid  and  hollow  portions  of  the  drive  shaft.  This  was 
accomplished  by  decoupling  the  airfoil  from  the  transmission  and  setting 
both  parts  of  the  system  at  their  maximum  noseup  positions  (e.g.,  to  set 
otM  =  lit  deg,  the  airfoil  was  placed  at  +  o’  =  lh  +  8  =  22  deg  and  the 
transmission  cam  was  adjusted  to  its  maximum  positive  displacement).  The 
airfoil  was  then  locked  to  the  transmission  by  a  clamp.  During  this  pro¬ 
cedure  a  second  angular  position  transducer,  mounted  on  the  solid  shaft, 
was  used  to  measure  the  input  angular  time  history. 

STEADY-STATE  TESTS 

The  steady-state  normal  force  and  moment  were  obtained  for  the  two  test 
airfoils  in  the  Acoustics  Research  Tunnel.  The  tunnel  speed  was  275  fps 
corresponding  to  a  Mach  number  of  0.25.  These  steady -state  data  were 
used  to  modify  the  tabulated  unsteady  data  of  Reference  1  to  obtain  the 
unsteady  data  for  this  study.  Furthermore,  the  Time  Delay  method  of 
analysis  (described  in  a  later  section)  uses  empirically  established 
variations  of  the  steady  data  to  model  unsteady  effects. 

The  unsteady  normal  force  and  moment  data  reported  in  Reference  1  were 
obtained  in  the  UARL  Low  Speed  Cascade  Tunnel,.  The  Reynolds  number  and 
Mach  number  of  the  current  test  program  were  approximately  the  same  as 
those  of  the  tests  of  Reference  1.  However,  the  Acoustics  Research  Tunnel 
used  for  this  study  is  an  open-jet  facility,  and  differences  were  expected 
in  the  steady  and  unsteady  aerodynamic  forces  and  moments. 
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The  steady  normal  force,  moment  and  drag  were  obtained  using  a  three-com¬ 
ponent  balance  and  standard  force  balance  techniques  typical  of  those  dis¬ 
cussed  in  Reference  2.  Figures  6  through  8  show  the  uncorrected  lift, 
drag  and  pitching  moment  versus  angle  of  attack  for  both  the  NACA  0012  air¬ 
foil  (top  panels)  and  SC  1095  airfoil  (bottom  panels)  as  obtained  in  the 
open-jet  wind  tunnel. 

In  Reference  2  Pope  discusses  the  effect  of  making  these  measurements  in  an 
open-jet  (i.e.,  no  floor  or  ceiling)  wind  tunnel  and  derives  the  following 
corrections  whereby  measurements  made  in  an  open- jet  can  be  modified  to  be 
consistent,  with  free  air  measurements: 


where  Aa,  Ac<q,  and  Ac-^^  are  the  amounts  by  which  the  original  measure¬ 
ments  must  be  modified  at  equivalent  values  of  c^_ ,  c  is  the  airfoil  chord, 
and  h  is  the  Jet  neight.  If  these  modifications  are  applied  to  the  data  of 
Figures  6  through  8  the  corrected  cq,  cq  and  c  versus  angle  of  attack 
curves  of  Figures  9  through  11  are  obtained.  As  a  check  on  the  accuracy  of 
these  results,  Figure  12  shows  the  comparison  of  the  corrected  drag  polar 
(i.e.,  c-l  vs.  cq)  data  for  the  NACA  0012  airfoil  to  the  results  on  page  9h 
of  Reference  3.  The  agreement  is  considered  to  be  good. 

DYNAMIC  TESTS 

After  the  steady-state  balance  tests  were  completed,  the  balance  was  re¬ 
moved  and  the  airfoil  was  attached  to  ti.e  oscillating  mechanism  previously 
described.  The  wind  tunnel  velocity  was  established  from  the  pressure  dif¬ 
ference  between  the  total  pressure  upstream  of  the  inlet  nozzle  and  the 
static  pressure  at  the  jet  exit.  The  wind  velocity  for  all  tests  was 
275  ±  2  fps. 

The  time-varying  input  angular  motion  from  the  transmission  and  the  output 
angular  motion  at  the  free  end  of  the  airfoil  were  obtained  from  the  angu¬ 
lar  position  transducers  and  were  recorded  on  FM  tape.  The  airfoil  and 
drive  mechanism  were  coupled  as  described  earlier.  The  airfoil  mean  angle 
of  attack  was  measured  with  a  clinometer.  Calibration  records  were  taken 
on  FM  tape  of  the  steady  zero  output  (airfoil  at  zxjvj),  the  steady  full-scale 
outputs  (airfoil  at  ±  8  deg)  and  the  unsteady  no  flow  time  history  at 
=  8  cps.  The  tunnel  velocity  was  then  set  to  275  fps  and  unsteady  time 
histories  were  recorded  at  nominal  driving  frequencies  of  8.0,  10.0  and 
12.5  cps.  This  gives  three  torsional  natural  frequency  ratios  (ug/fi)  for 
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the  same  torsional  natural  frequency.  This  test  procedure  was  repeated  at 
each  angle  of  attack  for  all  the  various  combinations  of  pitch  axis,  tor¬ 
sional  inertia,  torsional  natural  frequency  and  airfcil  type.  The  total 
matrix  of  test  conditions  is  shown  in  Table  I. 

ne  test  program  outlined  above  generated  such  a  large  volume  of  time- 
varying  data  that  cycle-by- cycle  inspection  and  data  reduction  were  not 
practical.  Accordingly,  the  data  reduction  was  automated,  with  the  bulk 
of  the  processing  being  done  on  the  Univac  1108  digital  computer.  The  FM 
tape  was  converted  to  a  digital  format  using  the  UARL  WISARD  data  system 
described  in  Reference  k  and  procedures  discussed  in  Reference  1.  Data 
were  recorded  for  sufficient  lengths  of  time  such  that  more  than  ten  cycles 
of  the  IP  driving  motion  were  included  at  the  lowest  oscillatory  frequency 
(8.0  cps ) .  A  typical  _,et  of  time  histories  for  the  basic  reference  condi¬ 
tion  (NACA  0012  airfoil,  25  percent  pivot  axis,  blade  alone  (reference) 
inertia,  and  5P  natural  frequency  ratio)  is  shown  for  the  four  mean  angles 
of  attack  in  Figure  13.  These  time  histories  were  obtained  directly  from 
the  computer  by  a  10-cycle  signal  averaging  process.  This  procedure  en¬ 
hanced  the  repeatable  elements  of  the  time  history  while  at  the  same  time 
reducing  the  noise  or  random  portion  of  the  signal.  The  signal  to  noise 
ratio  was  improved  by  ^/lO  =  3.16  by  this  process  (see  pages  Il6-ll8  of 
Reference  5). 

The  elastic  torsional  deflection  of  the  airf_.il  (the  difference  between 
the  total  angular  motion  and  the  input  angular  motion)  was  obtained  for 
each  test  condition  by  subtracting  the  input  angular  position  time 
history  from  the  averaged  airfoil  angular  position  time  history. 

0(t)  =  a(  x )  -|aM  +  .j-sin  2tttJ  (41 

where  0(x)  is  the  difference  between  the  measured  nondimensional  angular 
time  history  response,  ot(x),  and  the  input  driving  motion.  The  non- 
dimensional  time,  t,  is  given  by  t/P,  where  the  period,  P,  was  established 
from  the  ten-cycle  time-averaging  process  for  that  run.  Figure  l4  shows 
the  elastic  deflection  time  histories  corresponding  to  the  angle  of  attack 
time  histories  of  Figure  13. 

Two  measures  of  stall  response  are  noted  in  Figure  l4 .  The  first  para¬ 
meter  A0^  is  half  the  amplitude  of  the  first  cycle  of  unsteady  motion  which 
oo-responds  to  the  initial  stall.  The  second,  Oi^ppp,is  half  the  peak -to -peal 
or  maximum  deviation  angle.  These  two  parameters  were  extracted  from  each 
time  history  and  are  presented  in  Table  II  for  the  10  cps  driving  frequency. 
Since  stall  flutter  is  not  always  present  or  clearly  defined,  these  two 
parameters  could  not  always  be  extracted  and  are  not  entered  in  the  table 
for  such  cases. 

Examination  of  the  angle  of  attack  time  histories  indicated  that  there  were 
some  statistical  variations  in  measured  response  when  stall  flutter  was 
present.  It  was  recognized  that  the  signal  averaging  process  could  result 
in  significantly  reduced  unsteady  amplitude  if  there  was  appreciable  phase 
shift  from  cycle  to  cycle.  Sample  comparison  of  measured  cycles  of  elas¬ 
tic  deflection  angle  data  and  the  corresponding  time- averaged  cycles  are 
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shown  in  Figure  15.  The  top  panel  of  Figure  15  which  is  typical  of  most  of 
the  data  shows  that  little  statistical  variation  exists  and  that  the  time- 
averaged  cycle  is  representative  of  the  individual  cycles  of  motion.  For 
an  extreme  stall  flutter  condition  shown  in  the  lower  panel,  the  measured 
data  shows  a  mc^e  transient  characteristic,  and  unfortunately,  time¬ 
averaging  the  data  distorts  the  measured  response  during  the  period  of 
maximum  flutter.  The  general  character  of  the  response  (i.e.,  frequency, 
time  and  amplitude  of  initial  stall  response)  is  retained,  however.  Be¬ 
cause  this  phenomenon  was  not  typical  of  most  of  the  data  and  because  no 
obvious  way  of  choosing  the  average  or  nominal  cycle  of  data  for  each  of 
the  nearly  200  corditions  tested  was  evident,  the  time-averaged  cycle  of 
motion  was  used  in  the  time  history  and  amplitude  correlation  studies 
to  follow. 

There  are  many  ways  in  which  the  data  from  Table  II  can  be  presented  in 
graphic  form.  Figure  l6  depicts  the  initial  stall  angular  deviation  A  0^ 
versus  mean  angle  of  attack  for  the  NACA  0012  airfoil.  The  side  and  top 
headings  of  the  six  inserts  identify  the  particulai  insert  as  to  pitch 
axis  location  and  the  nominal  IP  frequency  2.  In  each  small  figure, 
there  are  four  curves  which  are  further  identified  to  indicate  the  moment 
of  inertia  and  the  nominal  airfoil  natural  frequency  (these  were  intended 
to  be  50  cps  and  70  cps  although  one  frequency  was  inadvertently  set  to 
55  cps).  Corresponding  results  for  the  SC  1095  airfoil  are  shown  in 
Figure  17.  Trends  obtained  in  these  tests  will  be  discussed  in  a  later 
section . 
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CORRELATION  STUDY  OF  TWO-DIMENSIONAL  RESULTS 


A  major  goal  of  the  two-dimensional  flexured  ainoil  test  program  was  to 
provide  dynamic  response  data  which  could  he  compared  with  predictions 
based  on  various  unsteady  aerodynamic  methods.  The  study  of  helicopter 
rotor  stall  flutter  experienced  in  flight  is  complicated  by  uncertainties 
as  to  rotor  inflow,  spanwise  flow,  and  blade  bending  and  plunging  motion. 
In  the  present  tests  these  effects  were  eliminated  and  direct  comparison 
of  unsteady  aerodynamic  models  was  possible.  Two  analytical  models  were 
considered  as  discussed  below. 

The  single  torsional  degree  of  freedom  differential  equation  of  motion  for 
the  flexibly  mounted  airfoil  section  oscillating  in  the  wind  tunnel  test 
section  is  given  by 

la  +  Ca  +  K(a  -  %)  =  M  (t)  +  T  sinftt  (5) 

11  a  o 


where  I 

airfoil  toisional  inertia  per  unit  span, 

(ft -lb-sec 2 /ft ) 

a 

= 

airfoil  angle  of  attack,  rad 

C 

= 

equivalent  mechanical  damping  per  unit  span, 

( ft-lb-sec/rad) /ft 

K 

= 

torsional  spring  constant,  (ft-lb/rad) /ft 

aM 

- 

mean  angle  of  the  oscillation,  rad 

Ma(t) 

= 

applied  aerodynamic  moment,  ft-lb/ft 

T 

0 

= 

magnit\;de  of  the  applied  torque,  ft-lb/ft 

O 

Ufa 

= 

angular  frequency  of  the  applied  torque,  rad/sec 

t 

= 

time,  sec 

The  numerical  techniques  used  to  solve  this  equation  are  detailed  in 
Appendix  III  of  Reference  1. 

DESCRIPTION  OF  THE  UNSTEADY  MODELS 

a ,  A,  B  Model 

In  the  present  investigation  two  models  of  the  unsteady  applied  moment, 
Ma(t),  were  examined.  In  the  a,  A,  B  method  described  in  Reference  1  un¬ 
steady  aerodynamic  moment  is  assumed  to  be  a  function  of  angle  of  attack 
and  its  first  two  time  derivatives.  Tabulations  of  aerodynamic  coefficients 

as  functions  of  a,  A  =  ca  and  B  =  (  ^Y  a  (where  c  is  the  airfoil 

2U0  \  2/  U^? 
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chord  and  UQ  is  the  free-stream  velocity)  were  generated  for  an  NACA  0012 
airfoil.  The  airfoil  was  mounted  in  a  wind  tunnel  and  oscillated  in  pitch 
at  combinations  of  mean  incidence  angle,  amplitude  of  angular  oscillation 
and  reduced  frequency.  Instantaneous  lift  and  pitching  moment  coefficients 
based  on  chordwise  pressure  distributions  and  instantaneous  values  of  a, 
a  and  o  were  determined  during  the  oscillating  airfoil  tests.  Curve 
fitting  procedures  were  then  used  to  determine  tables  of  lift  and  pitching 
moment  as  functions  of  o.  A,  B.  It  was  demonstrated  in  Reference  1  that 
unsteady  normal  force  and  moment  data  generated  during  the  sinusoidal 
airfoil  tests  could  be  used  to  predict  the  aerodynamic  response  of  an  air¬ 
foil  executing  a  nonsinusoidal  motion.  In  a  limited  number  of  flexured 
airfoil  tests  described  in  Reference  1,  good  qualitative  correlation  of 
measured  and  predicted  airfoil  dynamic  response  was  achieved.  To  model 
the  present  test  conditions,  two  changes  were  made  to  the  calculation 
procedure.  First,  sxnce  pitch  axis  was  a  variable  in  this  study,  pro¬ 
vision  was  made  to  include  pitching  moment  due  to  chordwise  offset  of  the 
aerodynamic  center  from  the  pitch  axis,  i.e.. 


CmxpA  A*  =  cmc/lt^a»  A»  B)  +  (*PA  “  0.25)03.(0,  A,  B)  (6) 

The  functional  terms  cmc/i|  and  c^  were  obtained  from  interpolation  of  the 
unsteady  data  tabulations  of  Reference  1. 


The  second  change  involved  variations  of  the  steady-state  lift  and  moment 
for  the  present  test  program  from  the  corresponding  steady  data  obtained 
in  Reference  1.  As  mentioned  earlier  in  this  report,  the  actual  steady 
force  and  moment  were  different  even  for  the  NACA  0012  airfoil  because  the 
tests  were  conducted  in  an  open  jet  wind  tunnel  and  because  the  airfoil 
effective  aspect  ratio  was  much  higher.  It  was  necessary  to  scale  the 
data  tabulations  of  Reference  1  to  account  for  these  effects.  The  method 
of  scaling  used  for  these  analyses  required  a  shift  in  the  entire  data 
tabulation  by  constant  values  of  angle  of  attack,  unsteady  lift  coefficient, 
and  moment  coefficient  according  to  the  following  relations. 


c-^a.  A,  B)  open  JET  =  Cl^a  +  6al*  A*  B^TAB  +  6cl  (?) 

and 

cm(a.  A,  B)  open  JET  =  °m^a  +  <5am»  A>  B-TAB  +  6cm  (8) 

The  constants  6a^,  60^,  6ci  and  5cm  were  established  for  each  airfoil  and 
were  equal  to  the  amount  of  shift  necessary  to  make  the  open  jet  steady- 
state  stall  points  (in  both  lift  and  moment)  match  the  steady-sxate  stall 
points  of  the  airfoil  of  Reference  1. 
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Time  Delay  Model 

Wine?  tunnel  airfoil  dynamic  response  was  also  calculated  using  the  recently 
developed  Sikorsky  Time  Delay  unsteady  aerodynamic  method.  This  formula¬ 
tion  has  been  developed  empirically  by  generalizing  the  results  of  a  set 
o:  oscillating  airfoil  test  progr.ims  and  is  aimeu  at  predicting  the  .1- 
steady  aerodynamic  characteristics  of  arbitrary  airfoils.  Its  goaf  ^s  to 
provide  the  blade  designer  with  unsteady  lift  and  pitching  moment  character¬ 
istics  of  various  airfoils  while  minimizing  the  need  to  conduct  extensive 
oscillating  airfoil  test  programs.  This  model  is  based  on  a  hypothesis  of 
physical  separation  process  and  differs  from  most  other  approaches  in  that 
it  does  not  depend  on  an  assumed  harmonic  variation  of  angle  of  attack. 

The  basic  assumption  is  that  there  exists  a  maximum  quasi-static  angle  of 
attack  at  which  the  pressure  distribution  and  the  boundary  layer  are  in 
equilibrium.  During  increases  in  angle  of  attack  beyond  this  static  stall 
angle,  there  are  finite  time  delays  before  a  redistribution  of  pressure 
causes  first  a  moment  breiik  and  then  a  lift  divergence  corresponding  to 
flow  separation.  The  relative  phasing  of  the  moment  and  lift  breaks  with 
angle  of  attack  produces  either  positive  or  negative  damping  of  the  motion. 

In  order  to  test  the  Time  Delay  hypothesis,  oscillating  airfoil  test  data 
from  Reference  6  were  examined.  It  was  noted  that  the  occurrence  of  dy¬ 
namic  lift  stall  could  not  be  correlated  with  maximum  angle  of  attack  or 
with  the  time  derivative  of  angle  of  attack.  The  onset  of  stall  was  found 
to  occur  prior  to,  coincident  with,  or  after  the  maximum  amplitude  of  the 
angle  of  attack  oscillation.  Best  correlation  of  the  angle  of  attack  for 
dynamic  lift  stall  was  achieved  by  evaluating  the  spread  between  the  static 
moment  stall  angle  and  dynamic  lift  stall  in  terms  of  elapsed  time  non- 
dimensionalized  by  free-stream  velocity  and  chord  length,  t*  =  AtgEp(U0/c ) . 
It  was  found  that  dynamic  lift  stall  generally  occurred  when  an  elapsed  time 
parameter  Tp  (which  starts  counting  when  the  static  stall  angle  is  reached) 
attains  a  t*  value  of  about  6. 


Dynamic  pitching  moment  stall  has  been  handled  similarly.  From  inspection 
of  the  test  data  it  was  evident  that,  in  general,  an  abrupt  change  in  the 
dynamic  moment  occurred  prior  to  the  lift  stall.  This  has  been  noted  in 
Reference  7  and  attributed  to  the  shedding  of  a  vortex  at  the  airfoil  lead¬ 
ing  edge  at  the  beginning  of  the  separation  process.  Rearward  movement  of 
the  vortex  over  the  surface  of  the  airfoil  tends  to  maintain  lift  but 
drastically  alter  the  pitching  moment.  Examination  of  pitching  moment 
hysteresis  loops  from  Reference  6  showed  that  pitching  moment  stall 
generally  occurred  when  the  elapsed  time  parameter  ip  reached  a  pitching 
moment  stall  time  constant  r_  of  about  2. 


In  order  to  apply  the  Time  Delay  model  to  a  given  airfoil,  only  static 
aerodynamic  data  are  required.  First,  the  airfoil  static  lift  and  pitching 
moment  data  are  used  to  define  the  approximate  variation  in  center  of 
pressure  between  the  static  moment  stall  angle  a  and  an  angle  of  attack 
above  which  the  center  of  pressure  is  assumed  fixed.  Second,  an  approxi¬ 
mation  is  made  to  the  c^  versus  a  curve  for  fully  separated  flow.  The 
sequence  of  events  occurring  during  one  stall-unstall  cycle  is  detailed  in 
Figure  l8.  Briefly  stated,  lift  and  pitching  moment  are  determined  from 


9 


potential  flow  theory  until  the  nondimensional  time  ip  reaches  x3  .  At  this 
point  the  pressure  distribution  begins  to  change,  leading  to  a  rearward 
movement  of  the  center  of  pressure  and  a  loss  in  potential  flow  pitching 
moment.  At  a  later  time  when  To  =  x#,  the  lift  breaks  from  the  static  line 
and  decreases  gradually  with  time  to  the  fully  separated  value,  CiSEP^' 
For  Tg  >  x*  the  center  of  pressure  coincides  with  C.P.gpp(a).  At  the 
point  where  a  =  0  the  rates  at  which  c^  approaches  cj  pfa)  and  C.P. 
approaches  C.P.ggp(a)  (if  it  does  not  already  equal  C.P.gEp(a))  are  in¬ 
creased.  When  a  falls  back  below  the  quasi-static  stall  angle,  a^,  the 
center  of  pressure  returns  to  the  quarter  chord,  potential  flow  pitching 
moment  effects  are  reintroduced  and  a  second  time  parameter  is  recorded 
to  govern  the  rate  at  which  c-^  returns  to  cqp^. 

Although  additional  correlation  studies  must  be  made  to  identify  the  effect 
of  airfoil  type  on  the  time  delay  constants  and  although  a  number  of  refine 
ments  to  the  present  model  may  be  implemented,  this  rather  simple  model 
gives  a  good  representation  of  the  essential  features  of  the  dynamic  stall 
process.  Correlation  typical  of  that  claimed  for  other  empirical  methods 
(References  9  and  10 )  has  been  found  with  data  from  References  1,  6,  and  8. 
Only  the  a,A,B  method  has  produced  better  correlation  (Reference  l),  but  it 
suffers  from  the  requirement  for  extensive  testing  and  data  processing. 
Figures  19  and  20  compare  the  MCA  0012  unsteady  lift  and  pitching  moment 
hysteresis  loops  measured  in  Reference  1  with  Time  Delay  results.  This 
correlation  was  achieved  by  setting  the  lift  break  time  constant  x*  equal 
to  k.O  instead  of  6.0.  Apparently, three-dimensional  effects  encountered 
in  this  test  reduced  the  time  interval  between  static  stall  and  dynamic 
lift  stall.  Also  shown  are  the  hysteresis  loops  predicted  using  the  a,A,B 
method.  The  a,A,B  method  correlation  is  with  the  data  from  which  the  a,A,B 
coefficients  were  derived.  In  addition  to  predicting  the  exact  form  of 
lift  and  moment  hysteresis  loops,  an  unsteady  model  should  faithfully  repre¬ 
sent  aerodynamic  pitch  damping.  Accordingly,  the  Time  Delay  model  was  used 
to  calculate  two-dimensional  aerodynamic  damping  for  the  reduced  frequency/ 
mean  angle  of  attack  test  points  of  Reference  11.  The  results  plotted 
versus  reduced  frequency  and  mean  incidence  angle  are  shown  in  Figures  21 
and  22.  Included  in  Figure  21  is  the  variation  in  damping  predicted  using 
potential  flow  theory.  Generally  excellent  correlation  of  measured  and 
predicted  damping  is  noted.  In  Figure  23, damping  calculated  with  the  Time 
Delay  model  and  the  a,A,B  method  using  the  data  tabulations  of  the  early 
unsteady  tests  (Reference  8)  are  compared.  The  accuracy  of  the  two  methods 
is  similar  when  compared  on  this  basis.  However,  as  noted  earlier  and  as 
shown  in  Figures  19  and  20  the  u,A,B  does  more  faithfully  reproduce  the 
original  hysteresis  loops. 

The  Time  Delay  model  was  applied  to  the  calculation  of  wi  nd  tunnel  airfoil 
dynamic  response  using  the  numerical  procedures  outlined  for  the  a,A,B 
solution.  The  uncorrected  steady-state  lift  and  pitching  moment  data 
shown  in  Figures  6  through  8  served  as  inputs  in  the  Time  Delay  calcula¬ 
tions.  Additionally  the  airfoil  mean  incidence  angle  used  in  the  Time 
Delay  solution  was  two  degrees  less  than  that  set  in  the  wind  tunnel.  The 
open  jet  flow  deflection  experienced  at  high  unsteady  lift  coefficients 
was  sufficient  to  decrease  actual  peak  angles  of  ttack  to  a  value  some¬ 
what  lower  than  the  geometrically  impressed  pitch  angle.  The  two  degree 
correction  to  gave  consistently  better  correlation  of  the  initial  stall 
time . 
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COMPARISON  OF  MEASURED  AND  PREDICTED  TWO-DIMENSIONAL  AIRFOIL  RESPONSE 


Correlation  between  measured  wind  tunnel  model  response  and  i  esponse  calcu¬ 
lated  with  the  unsteady  models  described  above  vas  examined  for  thirty-six 
test  conditions.  The  set  of  cases  studied,  denoted  in  Table  II,  was  suffi¬ 
cient  to  evaluate  the  independent  effects  on  airfoil  response  to  stall  of 
mean  incidence  angle,  torsional  natural  frequency,  chordwise  pitch  axis, 
torsional  inertia,  and  airfoil  type.  Relative  to  a  baseline  case  taken  to 
be  the  NACA  0012  airfoil  at  =  14°,  uiq/Q  =  5,  xp^  =  0.25,  and  I  =  Ippp, 
Figures  2b  through  27  show  measured  and  predicted  effects  of  mean  angle, 
torsional  natural  frequency ,  pitch  axis  and  airfoil  type  on  time  histories 
of  elastic  deflection.  In  Figure  24  the  measured  and  predicted  effects  of 
airfoil  mean  angle  of  attack  are  shown  for  the  NACA  0012  airfoil.  Compari¬ 
son  of  the  two  test  results  indicates  that  deeper  penetration  into  stall 
results  in  sharper  initial  stall  deflection  and  larger  residual  stall  flut¬ 
ter  response.  Each  unsteady  model  overpredicts  the  amplitude  of  response. 
As  shown  in  Figure  25  the  main  effects  of  an  increase  in  torsional  natural 
frequency  are  a  shift  in  response  frequency  and  a  decrease  in  the  ampli¬ 
tude  of  elastic  deflection.  The  lower  panel  of  Figure  25  shows  fair 
correlation  of  response  amplitude,  although  both  analyses  predict  an 
initial  stall  response  earlier  than  that  measured.  The  effects  of  moving 
the  airfoil  pitch  axis  forward  are  shown  in  Figure  26.  A  delay  in  the 
initial  stall  time  and  a  reduction  in  the  amplitude  of  response  are 
indicated  by  the  test  data.  The  analytical  results  do  predict  the  re¬ 
duction  in  response  amplitude,  but  the  Time  Delay  model  still  results  in 
overpredicted  response.  Finally,  Figure  27  compares  the  IIACA  0012  and  the 
SC  1095  airfoils  at  the  same  mean  angles  of  attack.  A  delay  in  the  initial 
stall  time  is  evident  for  the  SC  1095  airfoil  which  had  a  static  stall 
angle  measured  in  this  wind  tunnel  to  be  about  three  degrees  higher  than 
that  of  the  NACA  0012.  However,  the  SC  1095  stall  flutter  amplitude  was 
comparable  to  that  experienced  by  the  NACA  0012  at  this  condition. 

In  general,  the  time  history  correlation  was  acceptable.  The  point  and 
severity  of  initial  response  to  stall  and  the  response  frequency  were  well 
modeled.  The  trends  observed  in  test  were  matched  by  the  analysis, 
although  the  Time  Delay  model  generally  overpredicted  stall  flutter  re¬ 
sponse.  The  difficulties  encountered  in  time-averaging  the  stall  flutter 
response  data  may  have  contributed  in  some  instances  to  this  apparent  over¬ 
prediction  of  response.  The  basic  effects  of  structural  changes  on  blade 
response  time  histories  were  predicted  by  either  analysis. 

Wind  tunnel  and  analytic  results  were  also  compared  on  the  basis  of  initial 
stall  response  amplitude  and  overall  peak-to-peak  deflection.  The  time- 
averaged  data  were  reduced  and  it  was  found  that  the  initial  stall  response 
parameter  AO^  gave  the  most  consistent  indication  of  susceptibility  to 
stall  flutter.  The  possible  reduction  in  flutter  amplitude  introduced  by 
the  time-averaging  procedure  when  there  was  cycle-to-cycle  variation  in 
phase  made  it  somewhat  difficult  to  assess  the  amplitude  of  flutter  re¬ 
sponse.  Fortunately,  the  initial  stall  deflection  showed  virtually  no 
cycle-to-cycle  variation.  Tables  III  and  IV  show  the  response  amplitudes 
calculated  with  the  a ,  A,  B  and  Time  Delay  models  at  Q  =  10  cps .  Figure  28 
compares  the  measured  NACA  0012  deflection  angles  from  the  S3  =  10  cps  plots 
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of  Figure  l6  with  predicted  values  for  all  combinations  of  airfoil  torsion¬ 
al  natural  frequency  ratio,  torsional  inertia  and  pitch  axis.  The  corre¬ 
sponding  results  for  the  SC  1095  airfoil  are  shown  in  Figure  29.  Certain 
general  trends  of  deflection  angle  can  be  identified  in  both  the  test  and 
analytic  results: 

1.  Elastic  response  increases  with  mean  incidence  angle. 

2.  For  the  same  torsional  inertia,  the  response  is  generally 
greater  for  the  lower  frequency  airfoil  section. 

3.  The  amplitude  of  response  is  inversely  related  to  torsional 
inertia. 

4.  Forward  movement  of  the  pitch  axis  leads  to  a  decrease  in 
response . 

5.  The  SC  1095  airfoil  dynamic  stall  response  generally  begins 
to  build  up  at  a  higher  mean  incidence  angle  than  the  0012, 
but  the  two  airfoils  have  comparable  responses  once  stall  is 
penetrated. 

Although  the  amplitude  of  elastic  deflection  is  important  in  determining 
rotor  stability  and  performance,  the  torsional  moments  resulting  from 
stall  flutter  are  the  designer's  primary  concern.  To  gain  a  measure  of 
the  trends  of  torsional  moment  with  parameter  changes,  the  twisting  moment 
experienced  by  the  flexible  connector  in  the  model  airfoil  drive  system 
was  calculated  for  each  test  condition.  The  torsional  moment,  Mq,  was 
calculated  using  the  equivalent  spring  stiffness  of  the  connector. 

M0  =  Keq0  =  (I0  u)Q2)0  (9) 

where  Iq  represents  the  airfoil  torsional  moment  of  inertia.  The  torsion¬ 
al  moments  corresponding  to  the  AG.  values  of  Tables  III  and  IV  are  used 
in  Figures  30  through  33  to  show  tne  effects  of  blade  parameters  on 
structural  moments.  In  Figure  30  it  is  clear  that  decreasing  torsional 
natural  frequency  generally  reduced  stall  flutter  moments.  Although  the 
stiffer  system  experienced  lower  response  amplitudes,  the  corresponding 
structural  moments  (proportional  to  the  square  of  the  natural  frequency 
according  to  Equation  (9))  were  increased. 

As  seen  in  Figure  31  the  effect  of  forward  placement  of  the  airfoil  pitch 
axis  was  generally  to  decrease  vibratory  torsional  moments.  The  two  analy¬ 
ses  predicted  this  trend  with  comparable  accuracy.  That  forward  movement 
of  the  airfoil  pitch  axis  relative  to  the  aerodynamic  center  reduces  stall 
flutter  moments  can  be  understood  by  considering  lift  and  pitching  moment 
hysteresis  loops.  For  an  airfoil  with  pitch  axis  forward  of  the  center  o  1 
pressure,  positive  lift  forces  cause  negative  moments  about  the  pitch 
axis.  Usually  for  positive  lift,  the  lift  hysteresis  loop  is  traversed  in 
the  clockwise  direction,  which  results  in  a  pitching  moment  loop  (due  to 
lift)  which  is  traversed  in  the  counterclockwise  (stabilizing)  direction. 

A  decrease  in  test  torsional  moment  amplitude  with  decreasing  torsional 
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inertia  is  generally  seen  in  Figure  32.  The  a.  A,  B  program  shows  a 
somewhat  less  significant  trend.  Finally  the  two  airfoils  are  compared  in 
Figure  33.  The  trend  that  high  stall  flutter  moments  are  delayed  in  mean 
angle  with  the  SC  1095  airfoil  is  seen  for  two  different  combinations  of 
inertia  and  pitch  axis.  The  stalled  blade  loads  for  the  two  airfoils  are 
comparable.  Based  on  the  model  airfoil  test  results,  the  following  guide¬ 
lines  for  reducing  section  vibratory  torsional  moments  can  be  established. 

1.  Reducing  torsional  frequency  decreases  vibratory  moments 
although  increasing  the  amplitude  of  torsional  deflection. 

2.  Reducing  airfoil  section  torsional  inertia  increases  tin1 
elastic  response  amplitude  but  generally  decreases  staj.l 
flutter  moments. 

3.  Movement  of  the  airfoil  pitch  axis  forward  of  the  quarter  chord 
reduces  vibratory  moments. 

1*.  Use  of  the  cambered  airfoil  delays  the  onset  of  stall  flutter 
loads  to  higher  angles  of  attack  than  measured  with  the 
NACA  0012.  When  stall  is  penetrated,  however,  no  significant 
reduction  in  vibratory  moments  can  be  anticipated  with  the 
cambered  airfoil. 

These  general  trends  will  be  investigated  in  the  full-scale  blade  design 
study. 
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SURVEY  OF  EXISTING  STALL,  BOUNDARY  PREDICTION  TECHNIQUES 


As  an  introduction  to  the  study  of  helicopter  stall  boundaries,  a  literature 
survey  of  available  stall  boundary  prediction  techniques  has  been  conducted. 
The  objective  of  this  effort  was  to  evaluate  the  state-of-the-art  methodol¬ 
ogy  used  for  rapid  prediction  of  main  rotor,  stall-related,  operating 
boundaries  and  to  determine  whether  any  of  the  procedures  are  sufficiently 
comprehensive  to  be  of  general  usefulness  in  the  preliminary  design  process. 
Although  significant  progress  has  been  made  in  the  prediction  of  rotor  and 
control  system  loads  using  sophisticated  aeroelastic  analyses,  a  fairly  re¬ 
liable  but  less  time-consuming  approach  would  be  quite  valuable  in  prelim¬ 
inary  stages  of  design.  An  objective  of  the  present  study  was  to  determine 
whether  the  results  of  a  systematically  applied  design  study  conducted  with 
the  latest  aerodynamic  and  dynamic  modeling  techniques  can  be  used  to 
(l)  refine  the  simple  stall  boundary  prediction  techniques  or  (2)  offer 
insight  into  the  importance  of  blade  parameters  not  treated  by  current 
methods . 

Historically, it  appears  that  rotor  solidity  requirements  have  been  deter¬ 
mined  from  the  maximum  thrust  coefficient-solidity  ratios  required  in  hover 
and  forward  flight.  Forward  flight  has  been  the  more  difficult  part  of  the 
problem  to  analyze  because  of  the  complexity  of  the  operating  environment 
and  the  multiplicity  of  ways  in  which  stall  effects  can  be  manifested; 
e.g.,  buildup  in  power,  stresses,  gust  sensitivity,  vibration  and  control 
loads  or  loss  in  control  power  or  lifting  capability.  To  simplify  the  pre¬ 
liminary  design  process,  some  type  of  static  aerodynamic  stall  criteria 
has  generally  been  used.  The  aerodynamic  parameters  and  their  limiting 
values  at  the  onset  of  stall  have  been  chosen  through  extensive  correlation 
with  flight  test  data.  Obviously,  the  reliability  of  this  type  of  predic¬ 
tion  technique  decreases  with  the  extent  to  which  new  rotor  configurations, 
materials  and/or  flight  regimes  are  being  considered.  This  approach  also 
hampers  a  detailed  comparison  of  stall  boundary  prediction  techniques  found 
in  the  literature  because  individual  companies  have  generally  chosen  stall 
criteria  based  only  on  their  immediate  experience.  Detailed  comparison  of 
the  various  prediction  techniques  would  require  that  they  all  be  applied 
to  a  given  set  of  rotors.  This  comparison  would  require  more  specific 
information  about  the  manifestations  of  stall  than  is  generally  available 
in  the  literature. 

The  most  commonly  used  form  of  stall  boundary  is  a  plot,  of  maximum  thrust 
coefficient-solidity  ratio  versus  advance  ratio.  Figure  3^t  shows  some 
sample  stall  boundary  predictions  expressed  in  this  form  for  typical  pure 
helicopter  conditions.  Curve  1  derives  from  the  very  earliest  efforts  to 
predict  stall  boundaries  (Reference  12).  In  this  method  the  rotor  angle 
of  attack  at  the  tip  of  the  blade  at  270  deg  azimuth  is  calculated  using 
classical  rotor  theory  for  each  of  a  series  of  advance  ratios,  thrust 
coefficients  and  parasite  drag  coefficients.  Calculated  tip  angles  of 
attack  of  12  deg  and  l6  deg  are  taken  as  the  incipient  stall  boundary  and 
the  practical  operating  boundary.  Curve  1  taken  from  Reference  13  ap¬ 
parently  represents  the  incipient  stall  boundary  calculated  in  this  manner 
for  a  Bell  AH-1G.  Curve  2  is  a  stall  boundary  calculated  for  a  Kaman 
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UH-2C.  The  criterion  that  a  local  retreating  blade  angle  of  attack  exceed 
13.5  deg  was  stated  in  Reference  l4  to  agree  well  with  Kenan  experience. 
Blade  angles  of  attack  used  to  generate  this  curve  were  calculated  using 
the  NASA  charts  of  Reference  15 • 

Use  of  the  retreating  blade  angle  of  attack  as  the  criterion  for  the  rotor 
stall  limit  was  replaced  at  Sikorsky  by  a  method  based  on  the  correspondence 
of  the  knee  of  the  control  load  curve  with  a  limit  value  of  the  blade  pro¬ 
file  drag  torque  parameter  bCg  p/o  =  0.004.  Extensive  examination  of 
flight  test  data  for  different ’Sikorsky  type  rotors  showed  that  bCQ,o/a  ” 
0.004  indicated  incipient  stall  and  bCQ,p/a  =  0.008  corresponded  approxi¬ 
mately  to  the  practical  limit  of  operation.  Based  on  this  criterion,  stall 
boundaries  are  estimated  from  a  rigid  blade  trim  analysis  which  has  shown 
good  correlation  between  measured  and  predicted  rotor  performance.  In 
this  performance  program  the  effects  of  aerodynamic  design  factors  (e.g., 
number  of  blades,  chord,  airfoil  and  twist)  on  stall  boundaries  are 
modeled,  but  the  effects  of  blade  structural  design  properties  are  not  in¬ 
cluded.  Curve  3  shows  the  incipient  stall  boundary  calculated  for  the 
CH-53A  using  the  rotor  torque  criterion  bCQ,D/o  =  0.0C4. 

Boeing-Vertol  apnears  to  use  a  very  similar  approach  in  predicting  moment 
stall  inception.  Although  the  details  of  their  estimation  procedures  are 
not  available  ir  the  literature,  Reference  16  indicates  that  flight  test 
results  have  be<  n  generalized  and  a  limiting  value  of  a  moment  stall  param¬ 
eter  f(p,  ^R,  n,  v.-,  K,  up,  I@)  has  been  established  as  a  function  of  calcu¬ 
lated  angle  of  attack  at  270  deg  azimuth.  Based  on  this  limit  curve  the 
stall  boundaries  for  a  proposed  rotor  design  can  be  established  by  running 
a  rigid  blade  trim  analysis.  The  Boeing  procedure  appears  to  have  gone 
further  than  the  previous  methods  in  that  blade  torsional  frequency  and 
inertia  have  been  included  as  parameters. 

A  fourth  curve  included  in  Figure  34  is  taken  from  the  results  of  Reference 
17.  In  this  study  an  ideal  rotor  is  postulated  which  has  all  blade  sections 
operating  at  except  for  an  advancing  blade  region  for  which  c-^  is  re¬ 

duced  to  maintain  zero  rolling  moment.  This  calculation  is  based  on  an 
isolated  rotor  and  uses  steady-state  airfoil  data  and  uniform  inflow. 

The  results  shown  above  clearly  indicate  different  boundaries  for  different 
types  of  rotor  systems.  A  group  of  flight  test  boundaries  shown  in 
Figure  35  also  points  out  the  variety  of  stall  boundaries  presently  en¬ 
countered.  Curve  1  from  Reference  io  marks  a  stall  boundary  for  a  UH-1B 
rotor  and  is  apparently  determined  from  excessive  blade  torque  requirements 
and  oscillatory  chordwise  bending  loads.  Stall  flutter  is  not  a  problem 
under  these  conditions.  Reference  13  shows  the  overall  stall  boundary 
experienced  by  Bell  helicopters  which  is  characterized  by  rapidly  rising 
power,  blade  loads  and  controllability  problems.  Reference  10  indicates 
that  the  CH-47C  flight  envelope  is  limited  by  stall  flutter  pitch-link 
loads.  Curve  3  approximately  locates  this  stall  boundary.  The  stall 
operating  limit  measured  for  the  Kaman  101  rotor  is  also  shown;  however. 
Reference  19  dees  not  describe  the  form  of  stall  effect  encountered.  Fi¬ 
nally,  curve  5  depicts  the  boundary  corresponding  to  excessive  CH-53A 
vibratory  pushrod  loads . 
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In  summary,  the  results  of  the  literature  survey  indicate  a  variety  of 
stall  boundaries  and  preliminary  design  procedures.  None  are  documented 
in  detail.  However,  it  is  clear  that  structural  design  factors  as  well  as 
aerodynamic  factors  must  be  considered  if  generalized  stall  boundaries  for 
rotors  are  to  be  determined.  The  remainder  of  this  investigation  is  di¬ 
rected  toward  a  systematic  analytical  study  of  the  factors  contributing  to 
stall  boundaries  for  an  articulated  rotor.  This  study  is  preceded  by  a 
correlation  study  to  define  the  general  accuracy  of  the  aeroelastic  analy¬ 
sis  being  employed. 


FULL-SCALE  CORRELATION  STUDY 


A  full-scale  helicopter  correlation  study  was  performed  to  evaluate  the 
analytic  capability  for  predicting  the  effects  of  stall  on  rotor  character¬ 
istics  of  interest.  The  two  unsteady  aerodynamic  models  used  in  the 
correlation  of  two-dimensional  stall  flutter  results  were  also  compared 
based  on  full-scale  correlation.  It  was  essential  tl.  t  current  aeroelastic 
find  aerodynamic  modeling  techniques  be  evaluated  systematically  over  a  wide 
range  of  flight  conditions  so  that  confidence  in  the  aralysis  could  be 
gained  before  the  parametric  design  study  was  performed.  Additionally,  the 
correlation  study  served  to  point  out  problem  areas  which  require  improve¬ 
ments  to  the  current  analytic  methods. 

Control  loads,  blade  stresses  and  rotor  performance  were  calculated  for 
both  stalled  and  unstalled  flight  conditions  for  the  CH-53A  and  CH-5^B 
aircraft.  The  CH-53A  shown  in  Figure  3 6  is  a  high-speed  assault  transport 
helicopter.  It  has  a  rotor  diameter  of  72  feet,  a  chord  of  2.167  feet  and 
a  solidity  of  0.115.  The  nominal  tip  speed  is  696  ft/sec  and  the  -6  deg 
twist  blades  have  a  modified  NACA  0011  airfoil  section.  Nominal  aircraft 
gross  weight  is  3^,000  lb.  The  CH-5^B  is  shown  in  Figure  37.  The  CH-5^B 
rotor  is  essentially  the  same  as  that  of  the  CH-53A  but  with  twist  in¬ 
creased  to  -l6  deg.  Normal  gross  weight  is  38,000  lb.  Each  rotor  has 
6  blades.  Natural  frequencies  for  the  blades  of  the  two  aircraft  are 
shown  in  Table  V  for  the  tip  speeds  used  in  correlation  studies. 

ANALYTIC  METHODS 

The  Normal  Modes  Blade  Aeroelastic  Analysis  was  used  to  calculate  the  blade 
elastic  response  for  each  flight  condition.  In  this  analysis,  which  is 
described  in  Reference  20,  the  fully  coupled  blade  equations  of  motion  are 
solved  by  expanding  them  in  terms  of  uncoupled  flatwise,  edgewise  and  tor¬ 
sional  blade  modes.  The  modal  technique  facilitates  the  numerical  integra¬ 
tion  of  the  blade  equations  by  minimizing  dynamic  coupling  terms .  Given 
a  distribution  of  wake-induced  velocities  over  the  rotor  disc  and  an  aero¬ 
dynamic  model  to  represent  section  aerodynamic  loads,  blade  elastic  re¬ 
sponse,  blade  stresses  and  control  loads  are  calculated  for  each  fl'ght 
condition . 

Rotor  inflow  was  calculated  using  the  UAC  Prescribed  Wake  Aralysis  as  de¬ 
scribed  in  Reference  21.  The  function  of  this  program  is  to  compute  a 
rotor  circulation  distribution  that  is  compatible  with  a  prescribed  set  of 
blade  section  operating  conditions  and  a  prescribed  rotor  wake  geometry. 

In  the  present  efforts  a  helical  wake  structure  was  assumed.  Based  on  cal¬ 
culated  circulation  strengths,  blade  positions  and  wake  geometry,  this 
analysis  applies  the  Biot-Savart  Law  to  calculate  induced  velocity  at  each 
radial  and  azimuthal  blade  station.  Tc  avoid  unrealistically  high  values 
of  computed  velocity  when  blades  and  trailing  vortices  pass  close  to  one 
another,  velocity  inducing  effects  are  eliminated  for  blade-vortex  separa¬ 
tion  distances  less  than  a  specified  vortex  core  radius.  Because  the  blade 
response  program  employs  a  lifting  line  representation  of  blade  aerodynam¬ 
ics,  the  large  fluctuations  in  computed  inflow  velocities  which  result  from 
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assuming  too  smala  a  vortex  core  limit  lead  to  greatly  overpredicted 
fluctuations  in  lift  coefficients  and  blade  flatwise  bending.  In  this 
study  it  was  found  that  a  vortex  core  radius  of  one-half  the  blade  chord 
gave  good  agreement  of  predicted  and  measured  flatwise  blade  stress. 

An  iteration  procedure  is  used  to  couple  the  blad  '  response  and  circulation 
programs,  insuring  compatibility  of  the  induced  ve. ocities  and  the  blade  re¬ 
sponse.  For  a  given  flight  condition  the  normal  mode;'  program  is  first  run 
with  uniform  inflow  and  trimmed  automatically  to  specified  values  of  rotor 
lift  and  propulsive  force.  Tolerances  of  about  t 3 %  on  rotor  lift  and  ±8% 
on  propulsive  force  were  used.  One-per-rev  rotor  f looping  was  maintained 
at  less  than  2  -  j  degrees  by  adjusting  cyclic  pitch.  Based  on  required 
control  settings  and  resulting  motion,  rotor  induced  velocities  are  calcu¬ 
lated.  The  induced  velocity  distribution  is  applied  in  the  blade  response 
calculations  ana  the  rotor  is  retrimmed  by  adjusting  rotor  collective  pitch 
and  shaft  angle.  The  procedure  is  repeated  until  the  blade  motion  resulting 
from  the  response  calculations  is  consistent  with  that  for  which  the  in¬ 
duced  velocities  were  calculated.  Typi^ally^ two  iterations  are  required  to 
obtain  reasonably  converged  results. 

The  two  unsteady  aerodynamic  models  described  above,  namely  the  a,A,B  and 
the  Time  Delay  model,  were  used  to  calculate  blade  unsteady  lilt  and 
pitching  moment  coefficients.  In  the  absence  of  unsteady  drag  data  or 
proven  unsteady  drag  modeling  techniques,  steady-state  drag  data  were  used 
throughout.  A  smoothed  version  of  the  u,A,B  unsteady  data  tables  0“ 
Reference  8  was  used  in  the  correlation  study.  Results  of  Reference  1  in¬ 
dicated  that  this  set  of  unsteady  data  gave  best  full-scale  correlation 
results.  In  the  application  of  the  a,A,B  model  the  effects  of  blade  un¬ 
steady  motion  on  angle  of  attack  for  zero  lift,  lift  curve  slope  and  stall 
angle  were  applied  in  the  circulation  calculations. 

CH-53A  FLIGHT  TEST  CORRELATION  STUDY 

Extensive  flight  testing  of  the  CH-53A  aircraft  has  been  performed.  The 
primary  source  of  flight  test  data  chosen  for  the  correlation  study  is  the 
rotor  loads  program  described  in  Reference  22.  In  this  test  program,  de¬ 
tailed  measurements  of  blade  loads  and  responses  were  made  and  converted  to 
time  history  form.  Additional  vibratory  response  data  are  available  from 
the  results  of  Reference  23.  Based  on  the  data  of  Reference  23,  Figurp  38 
illustrates  the  buildup  of  vibratory  CH-53A  pushrod  loads  with  airspeed. 

A  rapid  growth  of  oscillatory  pushrod  load  is  noted  as  stall  is  penetrated. 
In  these  results  the  Sikorsky  ERITS  (Equivalent  Retreating  Indicated  Tip 
Speed)  procedure  has  been  used  to  correct  measured  airspeeds  for  small 
fluctuations  in  aircraft  gross  weight  and  load  factor,  tip  speed  and  air 
density.  ERITS  is  calculated  according  to  the  following  relation: 

E8ITS  « -  vc]yh>  u  do) 
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where  ftR  =  tip  speed,  kn 

Vq  =  calibrated  airspeed,  kn 
p  =  measured  density  of  air,  slugs/ft3 
pQ  =  reference  density,  slugs/ft3 
W  =  measured  gross  weight,  lb 
W0  =  reference  gross  weight,  lb 
LF  ~  load  factor,  nondimensional 

According  to  this  procedure,  described  in  Reference  2h ,  plotting  vibratory 
control  loads  versus  ERITS  reduces  the  scatter  in  experimental  da^a  and 
facilitates  comparisons  of  test  and  analytic  results.  When  the  mean  ERITS 
curve  has  been  compiled  for  a  given  aircraft,  mean  load  versus  airspeed 
curves  can  be  gt  erated  for  a  nominal  gross  weight,  tip  speed  and  density 
altitude. 

Correlation  of  control  system  loads,  blade  stresses  and  requii ed  power 
was  studied  at  a  nominal  aircraft  gross  weight  of  1+2,000  It  (Crj/a  =  0.083), 
a  tip  speed  of  710  ft/sec  and  a  3000  ft  density  altii  ude  f„r  airspeeds 
ranging  from  100  kn  to  170  kn.  The  effects  of  nonuniform  inflow  on  blade 

response  were  calculated  at  all  flight  conditions.  Also,  the  results 
obtained  with  the  two  unsteady  aerodynamic  models  were  compared  with  pre¬ 
dictions  using  steady  aerodynamics. 

The  inclusion  of  variable  inflow  was  found  to  be  essential  in  calculating 
the  proper  levels  of  blade  bending  moments.  Some  improvement  in  the  cor¬ 
relation  of  blade  torsional  moments  was  also  found.  Figure  39  shows  the 
blade  flatwise  and  edgewise  stress  correlation  obtained  with  a,A,B  unsteady 
aerodynamics  at  a  typical  stalled  flight  condition  (.137  kn/1+2,000  lb 
gross  weight).  Significant  improvement  in  the  correlation  of  blade  stress 
results  from  the  inclusion  of  variable  inflow  effects.  The  amplitude  of 
the  predicted  edgewise  stress  is  somewhat  less  than  the  measured  value.  It 
should  be  noted  that  the  CH-53A  blade  has  an  edgewise  natural  frequency 
near  1*P  which  undoubtedly  causes  computed  stress  to  be  sensitive  to  other¬ 
wise  small  errors  in  computed  1+F  loading.  Figure  1+0  compares  measured  and 
predicted  radial  distributions  of  vibratory  flatwise  and  edgewise  stress 
for  this  flight  condition.  The  improvement  in  stress  correlation  is  direct¬ 
ly  traceable  to  increased  third  quadrant  airloading  which  results  from 
higher  blade  angles  of  attack  calculated  by  the  variable  inflow  program. 

The  effect  of  inflow  on  the  correlation  of  pushrod  loads  is  indicated  in 
Figure  1*1.  Although  the  nonuniform  inflow  results  have  a  higher  predicted 
stall  flutter  amplitude^ the  triggering  of  the  oscillation  near  T  =  180  deg 
is  not  predicted.  Calculated  two-dimensional  angles  of  attack  are  general¬ 
ly  below  stall  at  this  position.  Possible  reasons  for  the  lack  of  agree¬ 
ment  in  the  stall  flutter  initiation  point  are  indicated  in  the  following 
section. 
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Because  it  is  fairly  well  established  that  accurate  prediction  of  stalled 
blade  loads  requires  the  treatment  of  unsteady  aerodynamic  effects, only  a 
few  comparisons  of  steady  and  unsteady  aerodynamics  were  made.  Figures  1+2 
and  1+3  compare  computed  blade  stresses  and  pushrod  loads.  The  increased 
edgewise  stress  calculated  with  steady  aerodynamics  apparently  reflects 
the  fact  that  without  dynamic  lift  effects,  higher  blade  angles  of  attack 
are  required  to  generate  a  given  rotor  thrust  resulting  in  higher  calcula¬ 
ted  drag  coefficients.  As  shown  in  Figure  1+3,  although  the  initial  calcu¬ 
lated  stall  moments  are  similar,  the  buildup  of  torsional  loads,  character¬ 
istic  of  an  undamped  or  negatively  damped  oscillation,  is  not  modeled  with 
steady  pitching  moment  data. 

The  a,A,B  and  Time  Delay  unsteady  aerodynamic  models  were  compared  at  the 
137  k.n/1+2,000  lb  flight  condition.  Although  Figure  1+1+  shows  comparable  com¬ 
puted  blade  stresses,  the  pushrod  loads  calculated  with  vhe  Time  Delay 
model  are  much  less  than  the  measured  values  (Figure  1+5).  Figure  1+6  shows 
similar  results  at  a  more  deeply  stalled  flight  condition  (155  kn/ 

1+2,000  lb  ).  Although  the  sharpness  of  the  initial  moment  stall  is  compar¬ 
able  with  the  two  unsteady  methods,  the  Time  Delay  results  do  not  indicate 
subsequent  oscillations  into  and  out  of  stall  with  pitching  moments  phased 
to  promote  the  growth  of  torsional  oscillations.  It  is  not  entirely  clear 
why,  relative  to  the  a,A,B  method, t;ie  Time  Delay  model  underpredicts  heli¬ 
copter  control  loads  while  overpredicting  the  stall  flutter  oscillations  of 
the  two-dimensional  wind  tunnel  model.  Examination  of  several  blade  sec¬ 
tion  pitching  moment/angle  of  attack  hysteresis  loops  indicates  not  so  much 
that  more  negative  pitch  damping  is  present  in  the  a,A,B  results  but  that 
pitching  moments  along  the  blade  are  more  in  phase,  leading  to  larger  modal 
excitation.  In  the  a,A,B  formulation, pitching  moment  coefficients  are  tab¬ 
ulated  as  functions  of  a, a  and  a.  When  a  torsional  response  is  present, 
blade  elastic  deflection  contributes  an  oscillation  to  the  ot,d  and  a  values 
all  along  the  blade  which  leads  to  more  similarly  phased  pitching  moments 
along  the  span.  In  the  Time  Delay  model  described  above,  moments  are  cal¬ 
culated  based  only  on  the  angles  of  attack  exceeding  the  steady-state  stall 
angle  for  a  certain  interval  of  time  and  are  independent  of  the  instantan¬ 
eous  angle  of  attack  time  histories.  For  small  differences  in  calculated 
angles  of  attack,  computed  pitching  moments  for  adjacent  blade  sections  can 
be  very  different  in  phase.  Because  the  wind  tunnel  airfoil  was  substan¬ 
tially  two-dimensional  in  nature,  the  effects  of  simultaneous  spanwise 
stalling  were  not  a  factor. 

That  better  correlation  of  stall  flutter  moments  was  possible  with  the 
a,A,B  method  is  evident  in  Figure  1+7,  which  shows  the  buildup  of  vibratory 
pushrod  load  amplitude  with  airspeed.  The  a,A,B  model  predicts  a  buildup 
rate  almost  identical  to  the  mean  of  the  test  data.  A  discrepancy  of  no 
more  than  10  knots  in  the  so-called  knee  of  the  control  load  curve  is  evi¬ 
dent.  Corresponding  results  for  flatwise  and  edgewise  vibratory  stress 
are  shown  in  Figure  1+8.  Maximum  vibratory  flatwise  stress  generally  oc¬ 
curred  at  about  the  0.8R  radial  station  and  maximum  edgewise  stress  oc¬ 
curred  at  0.5R.  Unfortunately,  blade  bending  stress  data  were  not  available 
at  the  highest  airspeeds.  Agreement  between  calculated  and  measured  vi¬ 
bratory  flatwise  stress  is  equally  good  with  the  two  models.  At  least 
until  deep  stall  is  reached,  the  analyses  (especially  the  Time  Delay 
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program)  generally  underpredict  vibratory  edgewise  stress.  Agreement  of 
measured  and  calculated  main  rotor  power  required  is  shown  in  Figure  1+9- 

CH-5to  FLIGHT  TEST  CORRELATION  STUDY 


The  stalled  flight  characteristics  of  the  CH-5^*B  rotor  system  were  studied 
recently  and  reported  in  Reference  2k.  Control  load  data  from  the  baseline 
flights  of  this  test  program  serve  as  the  basis  for  the  pushrod  load  cor¬ 
relation  study.  Blade  stress  and  required  power  data  were  taken  from  a 
general  CH-5^*B  flight  test  program  described  in  Reference  25.  The  CK-5^B 
flight  envelope  is  limited  by  the  buildup  of  vibratory  control  loads . 

Figure  50  illustrates  the  rapid  buildup  in  half-peak-to-peak  pushrod  loads 
experienced  at  an  aircraft  gross  weight  of  about  1+7>0G0  lb,  and  a  reduced 
rotor  speed  (96$  Np),  corresponding  to  a  nominal  thrust  coefficient-solid¬ 
ity  ratio  Crp/o  of  0.10.  Borne  uncertainty  exists  with  regard  to  the 'zero 
azimuth  reference  point  of  these  data.  Based  on  other  CH-5^B  pushrod 
load  data,  the  zero  azimuth  position  has  been  taken  as  the  last  nosedown 
peak  in  pushrod  load  before  relief  from  stall  causes  noseup  first  quadrant 
deflection.  Relative  to  the  CH-53A  data  a  decreased  stall  speed  and  a 
sharper  buildup  in  vibratory  load  with  airspeed  are  noted.  Also,  a  larger 
one-per-rev  component  is  evident  in  the  azimuthal  signature  of  the  pushrod 
load  resulting  apparently  from  the  high  blade  twist. 

Calculated  CH-5^B  control  loads  were  generally  less  than  measured  values . 
Figure  51  compares  pushrod  load  time  histories  at  a  flight  condition  of 
115  kn,  1*7,000  lb  gross  weight  and  9b%  normal  rotor  speed.  Relative  to 
comparable  CH-53A  results,  the  increase  in  IP  blade  pitching  momenL  is 
predicted,  but  the  amplitude  of  the  stall  flutter  oscillation  is  underpre¬ 
dicted.  In  order  to  probe  the  stall  boundary,  a  higher  speed  flight  condi¬ 
tion  (130  kn/l*7j000  lb)  was  simulated  analytically.  The  plot  of  vibratory 
pushrod  load  versus  airspeed  shown  in  the  top  panel  of  Figure  52  indicates 
that  a  definite  stall  boundary  is  predicted  by  the  analysis.  Relative  to 
the  CH-53A  calculations,  a  decreased  control  load  stall  speed  and  an  in¬ 
creased  rate  of  buildup  with  airspeed  are  clearly  predicted.  Again,  higher 
loads  are  computed  based  on  the  a,  T  B  model.  The  test  results  depicted 
in  Figure  50  indicate  that  the  rapid  buildup  of  peak-to-peak  loads  with 
airspeed  results  from  high  frequency  stall  flutter  oscillations.  The 
comparison  of  measured  and  predicted  stall  flutter  spike  amplitudes  shown 
in  the  lower  plot  of  Figure  52  indicates  that  the  a,  A,  B  results  reflect 
this  buildup  much  more  accurately  than  do  the  Time  Delay  calculations.  A 
comparison  of  calculated  vibratory  stress  with  test  data  from  Reference  25 
is  shown  in  Figure  5^*.  Measured  blade  stresses  are  well  below  critical 
levels  at  the  control  load  stall  speed.  Although  flatwise  stress  is  over¬ 
predicted,  edgewise  stress  agreement  is  very  good.  Compared  to  results 
obtained  for  the  CH-53A  a  greater  increase  in  vibratory  flatwise  stress 
with  increased  twist  is  predicted  than  was  measured. 

GENERAL  DISCUSSION  OF  CORRELATION 

In  general,  the  degree  of  correlation  presented  for  the  CH-53A  and  CH-5^*B 
was  considered  to  be  acceptable.  The  increases  in  blade  stress,  pushrod 
load  and  required  power  with  severity  of  stall  were  well  modeled.  Best 
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results  were  obtained  with  the  a.  A,  B  unsteady  aerodynamic  model  and  non- 
uniform  rotor  inflow. 

The  azimuthal  correlation  of  measured  and  computed  blade  flatwise  and  edge¬ 
wise  stress  was  much  better  than  that  of  the  pushrod  loads.  Although  ap¬ 
proximately  correct  amplitudes  of  vibratory  torsional  moments  were  pre¬ 
dicted,  the  analytic  results  do  not  generally  predict  torsional  oscillations 
building  up  at  f  =  180  deg.  Pushrod  loads  calculated  for  the  CH-53A  in 
Reference  1  indicated  the  potentially  strong  effect-  that  blade-vertex  inter¬ 
section  can  have  on  the  initiation  of  blade  torsional  oscillations.  As 
indicated  above,  a  large  vortex  core  diameter  was  employed  in  the  calcu¬ 
lations  of  induced  velocities  to  avoid  the  uncertain  effects  of  blade- 
vortex  proximity  on  calculated  flatwise  stress.  It  is  possible  that 
smoothing  the  rotcr  inflow  in  this  way  reduced  torsional  excitation.  Prob¬ 
lems  of  this  type  could  be  avoided  by  applying  a  lifting  surface  aerody¬ 
namic  model.  In  fact,  the  results  of  Reference  26  indicate  that  for  close 
blade-voi'tex  passage,  lifting  line  theory  badly  overpredicts  instantaneous 
blade  lift  while  underpredicting  pitching  moment.  Another  factor  which 
could  contribute  to  the  discrepancy  in  the  stall  flutter  initiation  point 
is  the  uncertainty  in  the  modeling  of  dynamic  stall  phenomena  at  moderate 
and  high  Mach  numbers.  The  present  unsteady  models  are  based  primarily 
on  generalizing  low  Mach  number  dynamic  stall  behavior.  Also,  in  view  of 
the  improvements  in  blade  stress  correlation  resulting  from  the  inclusion 
of  variable  inflow,  a  systematic  study  of  the  effects  of  wake  distortion 
on  blade  stresses  and  torsional  moments  should  be  performed.  The 
possibility  that  blade  torsional  oscillations  are  initiated  on  unstalled 
blades  through  motion  of  the  control  system  should  also  be  considered. 
Finally,  three-dimensional  flow  effects  and  unsteady  drag  modeling  may 
impact  the  initiation  of  blade  torsional  oscillations. 

The  failure  of  the  Time  Delay  blade  response  program  to  predict  full-scale 
blade  torsional  moments  was  disappointing.  The  correlation  of  hysteresis 
loops,  two-dimensional  aerodynamic  damping  and  wind  tunnel  airfoil  dy¬ 
namic  response  indicated  that  the  method  is  capable  of  predicting  stall 
flutter  response.  When  applied  to  the  three-dimensional  situation,  blade 
response  was  underpredicted,  apparently  because  simultaneous  stalling  of 
adjacent  blade  stations  was  not  promoted.  Work  is  presently  under  way  to 
evaluate  the  possibility  of  including  a  spanwise  stall  propagation  model 
in  the  representation  of  blade  aerodynamics.  Based  on  modeling  of  the 
flow  separation  process,  the  Time  Delay  method  can  easily  be  changed  to 
speed  up  local  separation  due  to  stalling  and  separation  of  adjacent  blade 
segments.  Some  evidence  of  this  type  of  behavior  can  be  found  in  blade 
pressure  coefficients  measured  in  flight,  but  quantitative  modeling  of 
stall  propagation  remains  to  be  established.  Because  the  need  for  a 
general  unsteady  model  not  dependent  on  extensive  oscillating  airfoil  tests 
is  obvious  and  because  Time  Delay  computing  times  and  storage  requirements 
are  much  less  than  required  by  the  a.  A,  B  model,  development  of  the  Time 
Delay  approach  will  continue.  For  the  present  investigation,  however, 
better  correlation  was  achieved  using  the  a.  A,  B  procedure;  it  was 
therefore  chosen  for  use  in  the  stall  boundary  generation  and  parametric 
design  studies  to  follow.  Because  test  data  are  available  over  a  wider 
range  of  airspeeds  and  gross  weights  ,  the  CH-53A  aircraft  was  selected  for 
use  in  the  stall  boundary  investigation. 


GENERATION  OF  STALL  BOUNDARIES  FOR  THE  CH-53A  HELICOPTER 


In  order  to  define  the  limitations  to  helicopter  flight  envelope  imposed  by- 
blade  stall,  the  manifestations  of  stall  were  examined  for  a  series  of  ad¬ 
vance  ratios  and  rotor  thrust  coefficients  along  the  stall  limit.  To  probe 
the  stalled  flight  behavior  of  the  CH-53A  rotor,  eight  flight  conditions 
were  simulated  using  the  Normal  Modes  Blade  Aeroelastic  Analysis  in  addi¬ 
tion  to  the  four  conditions  for  which  CH-53A  correlation  was  studied.  The 
a,A,B  unsteady  aerodynamic  model  and  variable  rotor  inflow  were  used 
throughout.  Relative  to  the  CH-53A  control  load  boundary  from  Figure  35  , 
Figure  5^  shows  the  flight  conditions  investigated.  At  each  of  three  Crp/o 
values  a  sufficient  range  of  advance  ratios  was  studied  to  show  the  effects 
of  stall  on  rotor  control  loads,  blade  stresses,  required  power  and  aircraft 
controllability.  Limiting  values  of  each  of  these  quantities  were  estab¬ 
lished  and  a  limit  speed  based  on  each  of  the  types  of  stall  losses  was  de¬ 
termined  for  each  Op/a.  Consistent  test  data  (i.e.,same  tip  speed,  fuse¬ 
lage  center  of  gravity,  ambient  air  density  and  temperature,  etc.)  were  not 
available  at  each  of  the  analytically  simulated  flight  conditions.  Agree¬ 
ment  with  available  experimental  data  was  examined  where  possible. 


The  nature  of  the  analytic  pushrod  load  buildup  with  advance  ratio  and 
blade  loading  is  shown  in  Figures  55  and  5 6.  Earlier  initiation  of  stall 
response  and  larger  vibratory  loads  are  noted  as  airspeed  or  rotor  thrust 
is  increased.  The  top  panel  of  Figure  57  shows  the  variation  of  half- 
peak-to-peak  pushrod  load  with  advance  ratio  for  the  three  thrust  coeffi¬ 
cient-solidity  ratios.  The  control  load  limit  criterion  was  taken  to  be 
a  vibratory  pushrod  load  of  ±2200  lb  which  corresponds  to  the  pushrod 
structural  damage  fatigue  limit  for  the  CH-53A.  Measured  and  calculated 
limit  speeds  based  on  this  criterion  are  compared  in  the  lower  portion  of 
Figure  57-  Actually,  the  control  load  boundary  is  characterized  by  high 
vibratory  loads  in  all  rotating  and  fixed  control  system  hardware.  For 
example,  the  loads  experienced  by  the  fixed  system  servomechanisms  build 
up  at  a  sharp  rate.  Servo  loads  were  calculated  in  the  normal  modes 
program  by  transferring  rotating  system  pushrod  loads  into  the  fixed 
system.  Figure  58  compares  measured  and  calculated  vibratory  loads  on  the 
fourth  quadrant  servomechanism  which  generally  experiences  the  most  severe 
load  buildup.  Calculated  fixed  system  loads  are  sensitive  to  the  degree 
of  convergence  of  the  rotating  pushrod  load  waveform.  A  higher  degree  of 
scatter  in  the  analytic  fixed  system  loads  is  therefore  present,  and  only 
an  approximate  faired  curve  is  shown  through  the  analytic  data  points. 

Good  correlation  of  fixed  system  loads  depends  upon  accurate  prediction  of 
not  only  the  amplitude  but  also  th_  harmonic  content  of  the  pushrod  loads 
and  is  therefore  more  difficult  to  achieve.  Some  discrepancies  are  noted 
in  the  stall  speed  for  the  lower  thrust  coefficients,  but  the  rate  of 
servo  load  buildup  is  well  modeled. 


In  order  to  determine  the  flight  limits  which  result  from  excessive  blade 
stress,  flatwise  and  edgewise  stresses  were  computed  for  each  flight  condi¬ 
tion.  Maximum  vibratory  edgewise  stress  was  consistently  found  to  occur 
at  0.5R  and  maximum  flatwise  stress  at  about  0.8R.  The  calculated  varia¬ 
tion  in  maximum  stress  is  shown  in  Figure  59-  As  discussed  ii.  the 
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correlation  section,  good  agreement  between  measured  and  calculated  blade 
stress  has  generally  been  found.  A  vibratory  stress  of  ±6 0C0  psi  has  been 
used  as  the  stress  limit.  Actually,  allowable  vibratory  stress  is  de¬ 
pendent  on  local  steady  stress  and  therefore,  varies  along  the  blade.  The 
±6000  psi  value  corresponds  roughly  to  the  infinite  life  limit  for  the 
aluminum  spar  under  typical  values  of  steady  stress  at  0.5R. 


Two  additional  manifestations  of  stall  are  shown  in  Figure  60.  The  first 
is  the  loss  in  available  control  moment  experienced  in  stall.  In  order  to 
assure  controllability  and  to  maintain  an  adequate  control  margin  for  the 
execution  of  maneuvers,  the  pilot  generally  maintains  the  capability  to 
apply  a  nosedown  pitching  moment  with  a  longitudinal  control  input.  In  ac¬ 
cordance  with  Military  Specification  MIL-H-8501A,  an  available  control 
moment  equal  to  ten  percent  of  the  control  moment  available  in  hover  has 
been  taken  as  the  limit.  For  each  of  the  computer  cases  described 
above  (and  shown  in  Figure  5^),  an  additional  case  was  run  in  which  the 
longitudinal  control  input  Bqg  was  changed  from  the  trim  value  to  the 
built-in  limit  of  the  aircraft.  In  each  case  the  change  in  ncsedown  pitch¬ 
ing  moment  was  noted  and  compared  with  the  calculated  head  moment  available 
in  hover.  The  top  panel  of  Figure  6 0  shows  the  variation  in  available 
moment  with  airspeed  for  the  three  gross  weights.  These  results  correspond 
approximately  to  neutral  center  of  gravity  operation;  although  as  noted, 
earlier  head  moments  were  not  trimmed  exactly  to  zero  for  each  of  the 
reference  conditions.  Examination  of  the  results  indicated  that  the  control 
power  derivative  APM/ABpg  remained  nearly  constant  even  out  to  200  knots. 

In  other  words, the  control  limitation  is  primarily  imposed  by  the  limit 
in  stick  travel  or  allowable  cyclic  pitch. 

The  last  stall  penalty,  the  buildup  of  required  main  rotor  power,  is  also 
shown  in  Figure  60.  Good  correlation  of  measured  and  computed  power  re¬ 
quired  was  shown  for  the  intermediate  gross  weight  in  Figure  1+9.  The  indi¬ 
cated  limit  value  of  6000  horsepower  corresponds  approximately  to  the  nor¬ 
mal  rated  main  rotor  power  available  from  the  CH-53A  engines  at  the  3000-ft 
density  altitude  being  used  in  the  calculations. 

The  stall  boundaries  calculated  above  are  collected  in  Figure  6l  where  max¬ 
imum  thrust  coefficient-solidity  ratios  based  on  each  of  the  stall  losses 
are  plotted  versus  advance  ratio.  The  calculated  boundaries  actually  fall 
very  close  to  one  another , being  separated  by  no  more  than  30  knots.  In 
genera]  agreement  with  flight  experience,  blade  stresses  are  first  to  reach 
critical  values  at  lower  Op/o  while  control  loads  tend  to  become  more  sig¬ 
nificant  at  gross  weights  and  altitudes  for  which  C>p/o  is  high.  The  con¬ 
trollability  bounlary  tends  to  be  of  least  concern.  Reference  to  Figure  57 
indicates  that  the  control  load  boundary  experienced  in  flight  occurs  at 
lower  airspeeds  than  have  been  calculated.  Accordingly,  extensive  flight 
data  detailing  stress,  power  and  controllability  boundaries  have  not  been 
available,  especially  at  high  C^/c.  Nevertheless,  the  good  prediction  of 
blade  response  phenomena  experienced  as  the  rotor  penetrates  stall  indi¬ 
cates  that  the  stall  mechanisms  are  being  reasonably  modeled  and  that  pre¬ 
dicted  changes  in  the  stall  boundaries  accompanying  blade  design  changes 
can  be  used  with  confidence. 
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Of  the  boundaries  considered,  the  control  load  limit  imposes  the  most  fun¬ 
damental  limitation  on  possible  expansions  to  the  flight  envelope.  Although 
Figure  59  indicates  that  blade  stress,  particularly  edgewise  stress,  shows 
a  significant  buildup  at  the  higher  speeds,  relatively  simple  design  changes 
can  be  employed  to  alleviate  stress  for  a  given  set  of  applied  bending 
moments.  For  example,  as  noted  earlier,  the  high  levels  of  CH-53A  edgewise 
stress  result  from  the  near  resonant  forcing  of  the  first  edgewise  mode. 
Flight  tests  of  modified  CH-53  blade  designs  have  shown  that  the  edgewise 
stress  problem  was  alleviated  by  detuning  the  edgewise  mode  from  the  UP 
resonance.  Also,  since  the  control  derivative  APM/AB^g  does  not  show 
severe  degradation  with  forward  speed,  additional  control  moment  could  ap¬ 
parently  be  obtained  by  increasing  allowable  cyclic  pitch.  Finally,  al¬ 
though  required  power  does  build  up  rapidly  with  blade  stall,  an  increase 
in  installed  power  would  permit  increased  forward  speeds.  The  flight 
envelope  determined  by  control  loads  is,  however,  more  difficult  to  expand. 
Beyond  the  knee  of  the  curve,  control  loads  build  up  at  such  a  rapid  rate 
that  strengthening  of  the  control  system  to  withstand  higher  loads  would 
result  in  only  small  increases  in  maximum  speed.  Changing  blade  design  to 
reduce  the  magnitude  of  the  applied  loading  is  one  fundamental  approach 
that  can  be  followed  in  an  attempt  to  reduce  the  control  load  buildup.  The 
potential  of  this  approach  is  examined  in  the  parametric  design  study  which 
follows . 


25 


PARAMETRIC  DESIGN  STUDY 


The  effects  of  blade  structural  design  changes  on  the  CH-53A  calculated 
stall  boundaries  were  examined.  Results  obtained  in  the  correlation  and 
boundary  generation  studies  indicated  that  the  mechanisms  by  which  pen¬ 
etration  into  stall  increased  rotor  loads  were  being  well  modeled  by  the 
analysis.  It  was  appropriate,  then,  to  investigate  for  a  typical  helicop¬ 
ter:  (l)  how  the  blade  response  mechanisms  in  stall  are  affected  by 

changes  in  blade  structural  design  parameters  and  (2)  whether  blade  design 
changes  can  be  made  which  delay  the  adverse  effects  of  stall  and  permit 
expansion  of  helicopter  flight  envelopes.  Particular  emphasis  was  placed 
on  the  control  load  boundary  which  often  presents  an  important  limitation 
to  possible  expansion  of  helicopter  operating  envelopes.  References  27 
and  2b  describe  two  specific  efforts  to  reduce  control  loads  by  changes 
in  blade  and  control  system  design.  In  the  present  work  a  more  general 
study  of  the  effects  of  blade  parameters  on  stall  boundaries  was  conducted. 
Examination  of  the  torsional  equation  of  motion  which  is  solved  in  the 
Normal  Modes  Aeroelastic  Blade  Response  Program  indicates  that  a  number  of 
parameters  can  have  fundamental  effects  on  blade  torsional  response. 
Equation  (89)  of  Reference  20,  which  describes  the  normal  modes  differ¬ 
ential  equations  of  motion,  shows  the  sources  of  blade  torsional  excitation. 
It  is  possible  to  identify  the  blade  torsional  structural  moment  (which  is 
directly  related  to  the  pushrod  load)  as  being  approximately  equal  to  the 
sum  of  an  applied  aerodynamic  moment,  a  torsional  inertia  moment,  a  moment 
due  to  the  product  of  shear  forces  and  blade  bending  deflections,  and  a 
moment  due  to  the  chordwise  offset  of  the  blade  center  of  gravity  from  the 
elastic  axis. 

Ms  =  Ma  +  Mj.  +  MSD  +  MCq  (II) 

A  group  of  blade  parameters  was  chosen  which  altered  each  of  these  sources 
of  torsional  moment.  Independent  changes  were  made  to  each  of  the  follow¬ 
ing  eight  parameters,  and  the  effects  of  each  of  the  terms  on  the  mechanism 
of  blade  torsional  response  in  stall  were  noted. 

(1)  Airfoil 

(2)  Torsional  mode  shape 

(3)  Blade  flatwise  and  edgewise  stiffness 

(1)  Chordwise  center  of  gravity  location 

(5)  Chordwise  elastic  axis  location 

(6)  Blade  twist 

(7)  Torsional  natural  frequency 

(8)  Torsional  inertia 
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In  order  to  concentrate  efforts  on  the  most  important  parameters,  an  ex¬ 
ploratory  study  was  first  conducted  at  a  previously  studied  flight  condi¬ 
tion  which  exhibited  large  amplitude  retreating  blade  torsional  oscillations. 
A  ll+O  kn/l*9»000  lb  gross  weight  flight  condition,  which  is  seen  in  Figures 
57)  59  and  60  to  be  near  the  control  load,  blade  stress  and  power  bounda¬ 
ries,  was  chosen  for  this  study.  Figure  62  shows  the  total  blade  root 
structural  moment  for  this  flight  condition  and  the  component  torsional 
moments  contributing  to  the  total.  Moments  are  nondimensionalized  accord¬ 
ing  to  the  conventions  of  Reference  20.  It  is  clear  that  blade  torsional 
moment  results  largely  from  the  combination  of  an  aerodynamic  pitching 
moment  and  the  resulting  torsional  inertia  moment.  The  inertia  moment  is 
almost  directly  out  of  phase  with  the  aerodynamic  moment  and  is  of  greater 
magnitude.  The  resulting  structural  moment  is,  therefore,  in  phase  with 
the  inertia  moment.  The  so-called  shear-deflection  moments  are  small.  The 
CH-53A  blade,  which  has  its  elastic  axis  and  center-of-gravity  at  approxi¬ 
mately  the  quarter  chord,  does  not  experience  moments  around  the  elastic 
axis  due  to  inertial  forces  acting  at  the  center  of  gravity. 

Independent  changes  were  made  to  each  of  the  eight  blade  parameters  ,  the 
blade  response  calculations  were  repeated  for  this  flight  condition,  and 
computed  values  «f  control  loads,  blade  stress  and  power  were  compared 
with  those  calculated  for  the  standard  blade.  In  general,  two  parameter 
values  were  studied,  one  above  and  one  below  the  corresponding  CH-53A 
value.  This  procedure  permitted  the  establishment  of  reliable  trends. 

In  each  of  the  modified  rotor  configurations,  the  rotor  inflow  distribution 
was  taken  from  the  corresponding  standard  blade  case  and  rotor  lift  and 
propulsive  force  were  trimmed  to  required  values  by  incrementing  collective 
pitch  angle  and  shaft  angle.  Based  on  the  results  of  these  preliminary 
calculations,  the  four  parameter  changes  wnich  had  the  most  beneficial 
effect  on  control  loads  were  chosen  for  further  study.  For  these  para¬ 
meters,  calculations  were  then  made  at  a  series  of  advance  ratios  and 
thrust  coefficients  to  define  modified  stall  operating  limits. 

AIRFOIL  PROPERTIES 

The  first  blade  property  to  be  studied  was  blade  airfoil  shape.  It  vas 
found  that  increasing  airfoil  static  stall  angle  had  a  clearly  beneficial 
effect  on  vibratory  control  loads.  Comparison  of  leading  and  trailing 
edge  stall  airfoils  did  not  show  a  significant  difference.  Optimization 
of  airfoil  contour  from  the  standpoint  of  dynamic  stall  behavior  is  a 
subject  deserving  a  great  deal  of  further  research.  Two-dimensional  dy¬ 
namic  stall  tests  of  the  type  described  above  and  fundamental  research 
into  the  effect  of  blade  contour  on  the  dynamic  stall  process  such  as 
described  ir.  Reference  28  must  be  done.  Attempts  to  investigate  fully 
'Me  effects  of  airfoil  shape  on  blade  behavior  in  stall  are  limited  by 
a  lack  of  unsteady  airfoil  data  or  a  proven  method  of  modeling  the  un¬ 
steady  characteristics  of  a  given  airfoil.  In  the  present  study  the 
effects  of  two  changes  in  airfoil  characteristics  on  stall  boundaries 
were  investigated  using  the  scaling  techniques  described  in  Reference  8. 
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Firs1-- ,  a  symmetric  airfoil  which,  relative  to  the  NACA  0012,  had  a  20  per¬ 
cent  increase  in  static  stall  angle  was  used  in  the  calculation  of  blade 
response  at  the  l40  kn/J+9»000  lb  flight  condition.  Unsteady  lift  and 
pitching  moment  data  for  this  hypothetical  airfoil  were  generated  by 
scaling  the  a.  A,  B  NACA  0012  data  so  that  the  static  stall  angle  and 
clMAX  werv  increased  by  20  percent.  Figure  63  compares  the  c  (a,  o,  o) 
and  cj_(a,  o,  o)  values  for  the  standard  and  scaled  airfoil  da?a.  Standard 
steady-state  drag  data  were  used.  The  pushrod  load  time  history  shows  that 
incrffsing  the  stall  angle  (as  might  be  expected)  gives  a  reduction  in 
peak-tc-peak  load  of  more  than  1+0  percent.  Reference  to  Figure  57  indi¬ 
cates  that  if  the  slope  of  the  pushrod  load  versus  airspeed  curve  were  not 
altered,  an  increased  allowable  airspeed  of  30-1+0  knots  would  be  realized 
with  this  airfoil.  Improvements  in  lift-to-drag  ratio  and  blade  stresses 
were  also  predicted. 

A  second  change  to  the  unsteady  aerodynamic  data  was  made  to  approximate 
the  pitching  moment  characteristics  of  a  leading  edge  stall  airfoil.  To 
accomplish  this  goal,  a  correction  function  6m(a)  was  added  to  the  NACA 
0012  data  which  increased  the  sharpness  of  the  pitching  moment  stall. 
Corrections  of  the  form  shown  in  Figure  61+  which  increased  positive 
pitching  moment  with  angle  of  attack  below  the  static  stall  angle  were 
added  to  the  NACA  0012  unsteady  data. 

CnL£(a,  A,  B)  =  cmj^CA  0012(a*  A»  B)  +  6m(a)  (l2) 

Airfoil  dynamic  stall  characteristics  as  reflected  by  the  unsteady  data 
tabulations  were  not  changed  but  the  snarpness  of  the  static  pitching 
moment  break  was  increased  as  shown  for  a  typical  Mach  number  in  Figure  6b. 
Standard  lift  and  drag  data  were  used.  Comparison  of  pushrod  load  time 
histories  shows  that  the  leading  edge  stall  airfoil  increases  the  overall 
amplitude  of  the  pushrod  load  but  does  not  have  a  strong  effect  on  re¬ 
treating  blade  stall  flutter  moments.  Changes  in  blade  stresses  and 
required  power  were  not  significant.  Because  unsteady  airfoil  data  in 
the  a.  A,  B  form  are  only  available  for  the  NACA  0012  airfoil  and  because 
of  the  uncertainties  involved  in  scaling  the  unsteady  NACA  0012  data  to 
represent  airfoils  with  different  stall  characteristics ,  the  effect  of 
airfoil  type  on  blade  response  was  not  studied  at  other  flight  conditions. 

TORSIONAL  MODE  SHAPE 

Independent  changes  to  blade  “orsional  mode  shape  which  were  also  con¬ 
sidered  had  only  secondary  effects  on  control  loads  at  the  flight  condition 
studied.  Changes  in  mode  shape  alter  the  radial  variation  of  elastic  twist 
angle  and,  therefore,  change  the  angle  of  attack  distribution  and  the  re¬ 
sulting  airloads  and  pitching  moments.  Because  torsional  natural  frequency 
and  torsional  inertia  were  assumed  to  be  fundamental  parameters,  independent 
changes  to  the  torsional  mode  shape  were  made  by  varying  blade  and  control 
system  stiffness.  For  the  standard  torsional  frequency  ratio  of  6. OP,  the 
two  limiting  values  of  torsional  mode  shape  (namely  those  corresponding  to 
rigid  blade/flexible  control  system  and  flexible  blade/rigid  control  system 
configuration)  were  calculated  using  the  natural  frequency  portion  of  the 
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normal  modes  program.  Substitution  of  the  two  limiting  mode  shapes  into 
the  blade  response  calculations  had  little  effect  on  blade  torsional  re¬ 
sponse.  Examination  of  the  two  cases  which  were  trimmed  to  the  same  rotor 
lift  and  propulsive  force  indicated  that  changing  the  radial  distribute -n 
of  elastic  twist  did  not  significantly  affect  blade  angles  of  attack  in  the 
area  of  the  three-quarter  radius,  where  the  largest  pitching  moment  loading 
is  experienced.  Blade  aerodynamic  pitching  moments  and  resulting  torsional 
structural  moments  were,  therefore,  essentially  unchanged.  At  least  for 
deeply  stalled  flight  conditions  torsional  mode  shape  does  not  have  a  sig¬ 
nificant  effect  on  control  loads. 

BENDING  STIFFNESS 

The  effects  of  blade  flatwise  and  e<  gewise  stiffness  on  control  loads  and 
blade  stresses  were  examined.  Changes  to  blade  bending  stiffness  had  vir¬ 
tually  no  effect  on  calculated  pushrod  loads.  Blade  bending  stresses  were, 
however,  reduced  in  some  cases.  Blade  bending  response  affects  section 
angles  of  attack  and  therefore  can  alter  blade  airloads  and  pitching 
moments.  Moreover,  a  torsional  moment  is  produced  by  the  blade  flatwise 
and  edgewise  shear  forces  which  act  on  the  deflected  blade.  It  is  shown  in 
Reference  20  that  the  importance  of  these  shear-deflection  effects  is  re¬ 
duced  as  local  blade  flatwise  and  edgewise  stiffnesses  approach  the  matched 
condition.  Three  stiffness  configurations  were  examined.  In  the  first  two 
cases,  blade  bending  stiffnesses  were  changed  by  factors  of  0.5  and  1.5. 

For  tne  reduced  stiffness  blade  the  first  flatwise  frequency  dropped  from 
2.72P  to  2.6UP  and  the  first  edgewise  frequency  dropped  from  3.91P  to 
3.25P.  At  the  lUO  kn/U9,000  lb  flight  condition  the  most  significant 
effect  on  blade  loads  was  a  four  to  one  reduction  in  the  magnitude  of  the 
UP  edgewise  bending  moment  at  the  critical  radial  station  caused  by  the 
separation  of  the  first  edgewise  frequency  from  the  UP  aerodynamic  loading. 
Overall  edgewise  bending  moment  was  decreased  only  by  about  20  percent  be¬ 
cause  the  3P  bending  moments  were  increased.  With  the  stiffer  blade  the 
first  flatwise  frequency  changed  from  2.72P  to  2.79P  and  the  edgewise 
frequer  cy  increased  to  U.U5P.  In  this  case  the  Uf  component  of  edgewise 
stress  was  reduced  by  20  percent  but  the  overall  stress  reduction  was  only 
about  four  percent.  In  neither  case  was  there  a  significant  effect  on 
vibratory  control  loads.  A  third  calculation  was  made  in  which  blade  flat¬ 
wise  and  edgewise  stiffnesses  were  matched  so  that  shear-deflection  mo¬ 
ments  were  eliminated.  This  change  had  only  a  small  effect  on  blade 
torsional  response.  This  result  was  expected  from  the  previous  calculation 
of  the  sources  of  torsional  moment  for  the  standard  blade  where  it  was 
shown  that  the  shear-reflection  moments  were  quite  small  at  this  flight 
condition  (Figure  62).  It  is  clear  that  for  this  blade,  torsional  response 
is  essentially  uncoupled  from  flatwise  and  edgewise  bending  and  that  the 
primary  source  of  torsional  excitation  is  the  aerodynamic  pitching  moment. 

ELASTIC,  CENTER-0 F-GRAVITY  AND  AERODYNAMIC  CENTER  AXES 

An  area  which  was  investigated  in  some  detail  was  that  of  blade  chordwise 
axis  positioning.  In  the  normal  modes  analysis  the  effects  of  displacing 
the  chordwise  center-of-gravity  and  aerodynamic  center  locations  from  the 
elastic  axis  are  modeled.  The  elastic  axis  and  pitch  axis  are  assumed  to 
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be  coincident.  Based  on  the  reference  blade  which  has  all  three  axes  at 
approximately  the  quarter  chord,  forward  and  aft  movement  of  the  center  of 
gravity,  aerodynjnic  center  and  elastic  axis  were  considered  independently. 
Of  the  cases  considered,  movement  of  the  coincident  elastic  and  center-of- 
gravity  axes  forward  of  the  blade  aerodynamic  center  had  the  most  benefi¬ 
cial  effect  on  vibratory  control  loads. 

Figure  65  shows  the  effect  on  pushrod  load  of  forward  and  aft  placement  of 
the  center- of-gravity  axis  relative  to  the  aerodynamic  center  and  elastic 
axes.  It  is  well  known  that  chordwise  offset  of  the  blade  center  of 
gravity  (c.g.)  from  the  elastic  axis  can  have  a  significant  effect  on 
blade  torsional  response.  For  example,  classical  flutter  and  static 
divergence  phenomena  depend  upon  the  c.g.  location.  Coupling  exists  be¬ 
tween  the  blade  flatwise  and  torsional  motions  such  that  for  a  blade  with 
its  c.g.  forward  of  the  elastic  axis,  rotor  centrifugal  and  inertial  forces 
produce  nosedown  (stabilizing)  moments  when  the  blade  moves  upward.  In  the 
first  analytic  case,  movement  of  the  center  of  gravity  aft  of  the  elastic 
axis  to  0.3c  led  to  a  subharmonic  (0.5P)  blade  response  involving  coupled 
flatwise  and  torsional  deflections.  Rotor  lift  oscillated  between 
1+7,000  lb  and  60,000  lb  and  there  were  substantial  increases  in  vibratory- 
head  moments  and  control  loads.  On  the  other  hand,  well-behaved  analytic 
responses  were  calculated  for  blades  with  0.265c  and  0.22c  center~of- 
gravity  locations.  A  gradual  trend  of  decreasing  vibratory  pushrod  load 
with  forward  movement  of  the  c.g.  from  0.25c  to  0.22c  is  noted  in 
Figure  65.  The  main  effect  of  moving  the  center  of  gravity  toward  the 
leading  edge  is  that  the  centrifugal  forces  introduce  a  steady  nosedown 
moment.  Neither  center- of-gravity  shift  had  a  major  effect  on  blade 
stresses  although  a  required  power  reduction  of  about  10  percent  was 
predicted  for  the  forward  c.g. 

The  second  axis  configuration  studied  involved  a  ±0.05c  shift  of  the 
elastic  axis  relative  to  the  coincident  aerodynamic  center  and  the  center 
of  gravity.  In  each  of  these  cases,  the  increment  to  the  aerodynamic 
moment  arising  from  the  offset  of  the  aerodynamic  center  and  elastic  axis 
was  approximately  balanced  by  a  change  in  the  inertial  moment  introduced 
by  c.g.  offset.  As  a  result,  the  high  frequency  control  loads  were  nearly 
unchanged. 

Offset  of  the  aerodynamic  center  from  the  coincident  elastic  axis  and  the 
c.g.  was  found  to  have  the  most  beneficial  effect  on  control  loads  at  this 
flight  condition.  The  two-dimensional  airfoil  results  described  earlier 
indicated  that  forward  movement  of  the  pitch  axis  (elastic  axis  for  blade) 
relative  to  the  aerodynamic  center  had  the  potential  of  decreasing  stall 
flutter  moments.  Although  the  effect  is  not  as  strong  for  the  rotor  as 
it  was  with  the  model  airfoil.  Figure  66  does  show  a  decrease  in  vibratory 
pushrod  load  of  about  20  percent  with  movement  of  the  elastic  axis  forward 
to  the  0.22c  position.  An  increase  in  steady  nosedown  moment  and  elastic 
deflection  angle  are  predicted  for  this  case. 

Because  there  war  a  reduction  in  torsional  moment  predicted  at  this  condi¬ 
tion  and  measured  on  the  wind  tunnel  airfoil  for  blades  with  pitch  axis 
shifted  forward  to  0.22c,  similar  calculations  were  performed  throughout 
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the  CH-53A  flight  envelope.  Figures  6 7  and  68  summarize  the  potential  ef¬ 
fect  of  the  forward  pitch  axis  on  rotor  loads.  Relative  to  the  standard 
blade  results,  modest  reductions  in  control  loads  and  required  power  and 
slightly  elevated  levels  of  flatwise  stress  are  evident.  Limiting  air¬ 
speeds  based  on  control  load,  power  and  stress  considerations  were  deter¬ 
mined  at  each  thrust  coefficient,  and  the  corresponding  stall  boundaries 
are  shown  in  Figure  69.  A  15-kn^f  increase  in  the  pushrod  load  limiting 
speed  is  predicted  at  the  highest  C^/o  .  However,  the  limiting  speed 
determined  by  flatwise  stress  is  reduced  by  about  15  knots  at  the  highest 

Qp/o. 


BLADE  TWIST 


A  similar  investigation  was  performed  to  determine  the  effects  of  blade 
twist  on  stall  boundaries.  Linear  twist  angles  of  -12  deg  and  0  deg  were 
used  at  the  reference  flight  condition.  Small  improvements  in  vibratory 
control  loads  and  required  power  were  calculated  for  the  high  twist  blades. 
Because  there  was  a  calculated  reduction  in  vibratory  load  and  because 
there  is  current  interest  in  high  twist  blcdes,  seven  flight  conditions 
were  simulated  using  the  -12 -deg  twist  blades.  Typical  time  history  com¬ 
parisons  are  shown  in  Figure  70.  At  the  137  kn/42,000  lb  gross  weight 
flight  condition  shown  in  the  top  insert  of  Figure  70,  the  increased  twist 
causes  a  nosedovr:  moment  at  ll+O  deg  azimuth.  This  nosedovn  moi. ent  is 
largely  caused  by  shear-deflection  effects.  The  increase  in  blade  twist 
makes  the  advancing  blade  tip  angle  of  attack  more  negative,  which  in¬ 
creases  the  tip  down  bending.  Drag  forces  acting  at  the  tip,  therefore, 
produce  a  nosedown  moment.  Because  the  rotor  is  only  lightly  stalled  this 
oscillation  does  not  grow.  In  the  fourth  quadrant,  the  increase  in 
built-in  twist  redistributes  blade  lift  and  decreases  blade  tip  angles  of 
attack  so  that  applied  aerodynamic  moments  and  oscillatory  response  are 
reduced.  At  a  higher  speed  flight  condition  (l80  kn/35»000  lb),  the 
second  quadrant  nosedown  moment  is  larger  and  serves  to  initiate  a  stall 
flutter  oscillation  which  grows  in  amplitude  to  such  an  extent  that  the 
potential  unloading  of  the  retreating  blade  tip  is  not  effective  in 
reducing  vibratory  pushrod  load.  Figure  71  illustrates  the  fact  that 
high  twist  generally  decreased  the  magnitude  of  the  pushrod  load.  At 
the  highest  airspeeds,  however,  the  effect  of  increased  twist  is  detri¬ 
mental.  At  all  thrust  coefficients  the  rate  of  control  load  buildup  with 
airspeed  is  higher  for  the  high  twist  blades.  Corresponding  results  for 
power  and  blade  stress  are  presented  in  figures  71  and  72.  Large  in¬ 
creases  in  flatwise  bending  stress  result  from  the  negative  tip  loading 
of  the  advancing  blade,  and  allowable  airspeeds  determined  by  flatwise 
stress  are  reduced  at  each  Cp/cr.  Six  additional  cases  were  analyzed  to 
determine  the  effect  of  twist  on  available  control  moment.  In  each  case 
the  longitudinal  cyclic  pitch  Bpg  was  set  equal  to  the  limiting  value  and 
the  available  nosedown  pitching  moment  was  noted.  The  results  show  that 
there  is  only  a  small  change  (Figure  72 )•  The  five  stall  boundaries  are 
shown  in  Figure  73.  Pushrod  loads  are  alleviated  up  to  an  advance  ratio 
of  about  0.45.  It  is  clear,  however,  that  flatwise  stress  becomes  the 
most  important  limitation. 
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TORSIONAL  FhaictOENCY 


The  model  airfoil  results  described  earlier  showed  that  reducing  torsional 
frequency  decreased  stall  flutter  moments.  Additionally,  experimental 
and  analytic  studx^s  described  in  Reference  27  show  the  significant  impact 
of  torsional  frequency  on  control  loads.  The  results  of  Reference  27 
showed  that  reducing  torsional  frequency  either  by  reducing  blade  or 
control  system  stiffness  or  by  increasing  torsional  inertia  reduced  stall 
flutter  moments  uniformly.  It  was  verified  in  the  torsional  mode  shape 
calculations  above  that  for  the  same  torsional  frequency,  the  relative 
values  of  blade  and  control  system  stiffness  did  not  have  a  strong  effect 
on  stall  flutter  response.  For  the  present  work,  changes  in  torsional 
frequency  were  accomplished  by  varying  blade  and  control  system  stiffness 
proportionately  so  that  mode  shape  was  not  altered.  Relative  to  the  stand¬ 
ard  blade  frequency  of  6. OP,  blade  frequencies  of  3.0,  4.0,  5.0,  8.0,  12.0 
and  00  times  rotor  frequency  were  analyzed  at  the  l40  kn/49,000  lb  flight 
condition.  Figure  74  compares  the  pushrod  load  time  histories  calculated 
with  four  values  of  torsional  frequency.  Reducing  frequency  from  the  6. OP 
value  reduces  vibratory  pushrod  load.  The  largest  ’  ibrating  pushrod  loads 
occurred  for  the  12P  blades.  However,  the  12. OP  results  may  be  subject  to 
some  erroi  because  over  a  limited  azimuth  range  the  nondimens ional  u  and  a 
values  associated  with  the  12. OP  response  exceed  the  limits  of  the  unsteady 
data  tabulations.  The  stall  flutter  spike  amplitude,  shown  in  Figure  74, 
and  the  overall  amplitude  of  pushrod  load  are  plotted  versus  torsional  fre¬ 
quency  ratio  in  the  top  inserts  of  Figure  75.  The  elastic  torsi  onai  tip 
deflections  corresponding  to  the  stall  flutter  spike  and  to  the  overall 
vibratory  load  are  shown  in  the  lower  panels  of  Figure  75.  It  is  interest¬ 
ing  to  note  the-  at  least  in  the  frequency  range  from  3.0  to  6. OP  the  tor¬ 
sional  deflection  amplitude  corresponding  to  the  stall  flutter  spike  re¬ 
mains  relatively  constant. 

The  stall  flutter  spiking  frequency  is  compared  with  the  torsional  natural 
frequency  in  the  top  panel  of  Figure  76.  The  response  frequency  is  in  all 
cases  greater  than  the  natural  frequency,  and  the  spread  between  the  two 
increases  as  natural  frequency  is  reduced.  The  shift  in  frequency  can  be 
attributed  to  an  aerodynamic  stiffening  caused  by  the  generally  negative 
slope  of  the  pitching  moment  versus  angle  of  attack  curver .  The  aerody¬ 
namic  pitching  moment  acts  to  apply  a  restoring  moment  in  parallel  with  the 
structural  restoring  moment.  (For  negative  slopes  o  the  cm  versus  a 
curves,  as  the  blade  undergoes  noseup  torsional  deflection,  angle  of  attack 
is  increased  and  a  nosedown  pitching  moment  is  experienced. )  The  lower 
panel  of  Figure  j6  illustrates  this  effect  schematically.  The  hysteresis 
1oop  determined  by  the  generalized  aerodynamic  pitching  moment  and  the 
first  torsional  mode  response  is  shown  over  the  range  of  blade  azimuth  cor¬ 
responding  to  stall.  As  shown  in  the  figure,  the  variation  of  general¬ 
ized  pitching  moment  with  generalized  torsional  modal  response  0^  can 

be  linearly  approximated  in  terms  of  constants  c^  and  C2  oy 

MAg  =  C1  +  C2°G  (13) 

where  the  constant  Cp  is  negative.  Based  on  this  approximation  a  stall 
flutter  spiking  frequency  can  be  calculated.  For  blades  Vvithout  c.g. 
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offset,  the  differential  equation  used  to  solve  for  the  first  torsional 
modal  response  can  be  approximated  by 

\  =  !T  %  (J 

w 


where  Oq  represents  the  second  derivative  of  0^  with  respect  to  non- 
dimensional  time  fit  and  1^  denotes  the  generalized  torsional  blade  inertia. 
If  Maq  can  be  represented  by  Equation  (13). 
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so  that 


0n  +  (o)Qz  -  c2/Ig)0g  =  71 


and  the  torsional  response  frequency  is  given  by 


oj  =  fa  2  -  co/I„ 

response  v  0  ^  G 


Examination  of  the  calculated  generalized  moment /torsional  response 
hysteresis  loops  shows  that  with  the  present  airfoil  data  tabulations  a 
value  of  C2/l0G  approximately  equal  to  -2h/rad  is  found.  This  value  in 
Equation  (l8)  gives  response  frequency  as  a  function  of  natural  freauency, 
as  shown  in  the  top  panel  of  Figure  76.  This  figure  indicates  that  the 
trend  exhibited  by  the  normal  modes  cases  is  reasonably  well  represented 
by  this  simple  model. 

In  order  to  understand  the  beneficial  effect  which  reducing  frequency  has 
on  vibratory  control  loads,  the  analytic  solution  to  the  torsional  response 
equation  was  examined  for  the  wq  =  h.OP  blade.  Figure  77  shows  the  toi- 
sional  moments  experienced  by  the  soft  blade  and  the  standard  6. OP  blade. 
The  amplitudes  of  the  aerodynamic  pitching  moments  are  approximately  equal. 
It  is  clear  that  the  reduction  in  structural  moment  results  from  a  decrease 
in  the  inertial  moment  Mj.  (For  a  second-order  torsional  system2having 
inertia  I  and  responding  harmonically  at  a  frequency  wq  ,  Mj  =  Iu)q0),  The 
inertial  moment  is  greater  than  and  approximately  ou>  of  phase  with  the 
aerodynamic  moment,  so  that  for  a  given  amplitude  of  M, .  reducing  Mj  reduces 
the  overall  structural  moment. 

Explanation  of  the  mechanism  by  which  reducing  frequency  reduces  structural 
moments  requires  understanding  of  the  interaction  between  high  frequency 
blade  torsional  motion  and  applied  aerodynamic  moments.  The  high  frequency 
component  of  the  unsteady  aerodynamic  moment  depends  on  the  amplitude  and 
frequency  of  the  blade  elastic  response.  The  elastic  response  depends  in 
turn  on  the  applied  moment.  In  the  present  computation  procedure,  the 
pitching  moment  characteristics  are  extracted  from  the  c  (a,  A,  B)  tables 
where  the  a,  A  and  B  values  are  strongly  dependent  on  elastic  response.  The 
interpolated  pitching  moments  are  incorporated  into  the  timewise  integration 
of  the  blade  equations  of  motion  which  proceeds  until  compatible  blade 
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torsional  response  and  applied  aerody  imic  moment  time  histories  are  es¬ 
tablished.  Detailed  analysis  of  blade  response  results  is  made  difficult 
by  the  complex  dependence  of  aerodynamic  moments  on  unsteady  motions  and  by 
the  fact  that  the  different  blade  designs  iterate  to  different  solutions  in¬ 
volving  different  regions  of  the  ot,  A,  B  tables.  The  results  of  the  fre¬ 
quency  study  indicate  that  over  the  range  of  frequencies  from  3. OP  to  6. OP 
the  shapes  of  the  pitching  moment  coefficient /angle  of  attack  hysteresis 
loops  for  outboard  blade  sections  are  very  similar.  In  other  words  the 
amplitudes  of  the  aerodynamic  forcing  functions  extracted  from  the  a,  A,  B 
tables  are  relatively  insensitive  to  frequency  over  this  range.  Instead, 
the  aerodynamic  moments  are  amplitude  dependent  and  a  fixed  amplitude  limit 
cycle  response  is  reached  independent  of  frequency.  For  a  given  amplitude 
deflection,  blade  stiffness  accounts  for  the  difference  in  structural  loads. 
The  softer  blade  experiences  a  reduction  in  retreating  blade  vibratory  load 
of  about  50  percent  at  this  flight  condition.  For  frequencies  above  about 
6. OP  the  form  of  the  c  (a.  A,  B)  variation  changes  and  different  hysteretic 
behavior  is  evident,  in  this  frequency  range  the  dependence  of  the  aero¬ 
dynamic  response  on  frequency  is  stronger  and  calculated  moments  are  phased 
to  reduce  the  amplitude  of  the  torsional  response.  The  vibratory  pushrod 
load  and  elastic  deflection  results  shown  in  Figure  75  show  that  relative 
to  the  standard  blade,  deflections  for  the  8.0  and  12. OP  blades  are  de¬ 
creased  but  increased  stiffness  causes  comparable  twisting  moments. 

The  uTr  =  1+.0P  blade  was  selected  for  investigation  throughout  the  flight 
envelope.  For  this  blade  the  GJ  and  root  spring  values  were  0.1+5  X  the 
standard  values.  Figures  78  and  79  summarize  the  analytic  results.  Siz¬ 
able  reductions  in  the  overall  vibratory  pushrod  load  are  shown  in 
Figure  78.  Reducing  frequency  has  the  most  beneficial  effects  at  high 
thrust  coefficients.  Results  at  the  l60-kn,  high  C^/o  flight  condition 
may  be  subject  to  question  because  at  this  condition  the  retreating  blade 
is  so  deeply  stalled  that  angles  of  attack  exceed  the  limits  of  the  un¬ 
steady  data  tables.  Allowable  airspeed  increases  on  the  order  of  15  knots 
are  shown  for  the  pushrod  loads  at  the  two  highest  Cp/a  values.  The 
reduction  shown  in  required  power  results  from  the  fact  that  the  increased 
compliance  of  the  torsionally  soft  blades  permits  the  blade  to  twist  out 
of  deep  stall  more  readily  and  achieve  a  more  optimum  angle  of  attack 
distribution.  Also,  because  the  soft  blades  generally  experienced  a 
larger  steady  nosedown  elastic  twist  angle  for  the  same  rotor  lift,  the 
center  of  lift  was  shifted  inboard  and  tip  angles  of  attack  and  drag 
forces  were  reduced.  The  stress  results  (Figure  79)  indicate  only  a 
small  impact  of  frequency.  A  limited  number  of  calculations  were  made  to 
determine  the  effects  of  torsional  frequency  on  the  available  control 
moment  boundary.  A  small  decrease  in  control  moment  stall  speed  is  shown 
in  Figure  79-  Some  loss  in  pitching  moment  capability  may  result  from 
inboard  movement  of  the  blade  lift  forces  due  to  the  change  in  angle  of 
attach  distribution.  The  changes  in  stall  boundaries  are  presented  in 
Figure  80.  The  pushrod  load  results  show  increases  in  allowable  speed 
or  allowable  payload  up  to  an  advance  ratio  of  about  0.1+7  with  the  log  =  1+  .OP 
blades . 
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TORSIONAL  INERTIA 


The  final  parameter  to  be  investigated  was  blade  torsional  inertia.  Be¬ 
cause  inertia  changes  independent  of  torsional  frequency  and  mode  shape 
were  desired,  blade  inertia  and  blade  stiffness  and  root  stiffness  were 
varied  in  proportion.  Physically,  this  study  would  treat  blades  of  the 
same  material  and  section  geometry  but  with  differing  cross  section  scale. 
Because  torsional  response  phenomena  were  the  primary  concern,  blade  mass 
and  flatwise  and  edgewise  inertias  were  not  changed.  Three  inertia  values 
were  used  to  calculate  response  at  the  llO  kn/l9»000  lb  flight  condition. 
The  resulting  pushrod  load  time  histories  are  shown  in  Figure  8l.  In  a 
manner  similar  to  that  found  for  the  low  frequency  blades,  decreasing  the 
inertia  reduced  the  inertia  moment  without  significantly  affecting  the 
aerodynamic  moment.  A  limit  cycle  solution  whose  amplitude  was  rela¬ 
tively  insensitive  to  inertia  was  calculated.  The  decrease  in  vibratory 
moment  with  decreasing  inertia  at  the  same  frequency  is  not  in  agreement 
with  the  results  of  Reference  27  where  it  was  stated  that  the  controlling 
parameter  was  torsional  frequency  and  that,  for  a  given  frequency,  pitch- 
link  loads  were  not  dependent  on  the  choice  of  stiffness  and  inertia 
values.  To  clarify  this  situation,  cases  were  run  in  which  frequency 
(and  mode  shape)  were  altered  by  changing  inertia  only.  The  results  of 
these  cases  shown  in  Figure  82  imply  that  at  this  flight  condition 
(l40  knA9>000  lb  gross  weight),  decreasing  frequency  by  increasing  inertia 
is  less  beneficial  than  increasing  frequency.  As  a  summary,  the  amplitudes 
of  the  vibratory  pushrod  loads  from  all  the  frequency  and  inertia  cases  are 
plotted  versus  torsional  natural  frequency  ratio  in  Figure  83.  Distinct 
trends  are  obtained  by  changing  GJ,  GJ  and  inertia,  and  inertia  only.  The 
results  are  also  plotted  against  the  stiffness  ratio  (GJ/GJgrpp).  It  is 
clear  that  reducing  stiffness  below  the  standard  CH-53A  value  reduces 
vibratory  loads  regardless  of  whether  the  frequency  is  reduced.  It  is 
also  clear  that  the  independent  effect  of  reducing  inertia  is  to  decrease 
the  amplitude  of  suall  flutter  spikes. 

To  determine  the  effect  of  inertia  on  control  loads  for  the  same  torsional 
frequency,  the  behavior  of  blades  having  a  50-percent  reduction  in  in¬ 
ertia  and  GJ  was  examined  along  the  stall  boundaries.  The  results  are 
summarized  in  Figures  8^t  through  86.  Decreases  in  peak-to-peak  pushrod 
loads  were  predicted  at  all  flight  conditions.  Relative  to  the  1.0P  blade 
which  also  had  GJ  -  ^GJsqip),  the  results  of  Figure  8U  indicate  that  re¬ 
ducing  inertia  by  50%  generally  reduced  loads  even  further.  Pushrod  load 
amplitude  was  reduced  to  such  an  extent  that  extrapolation  of  the  results 
to  the  limiting  load  level  was  difficult  without  repeating  calculations  at 
airspeeds  higher  than  previously  studied.  Because  blade  stress  and  power 
limits  had  already  been  reached,  this  was  not  done  and  the  limiting  air¬ 
speeds  are  determined  only  approximately.  Reductions  in  required  power 
comparable  to  those  calculated  with  the  Wq  =  1.0P  blades  were  predicted. 
Gmail  degradations  in  blade  stress  and  available  control  moment  are  shown 
in  Figure  85  for  the  low  inertia  blades.  The  stall  boundaries  based  on 
each  of  the  five  performance  indicators  are  shown  in  Figure  86. 

The  results  of  the  anaiytic  parameter  variation  study  can  be  used  to  pro¬ 
vide  the  following  set  of  general  design  guidelines  aimed  at  reducing  high 
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frequency  control  loads.  Findings  1-1+  apply  most  directly  to  low  twist 
blades  with  symm^tri ;  J.  airfoil  sections  and  having  little  aerodynamic 
pitching  moments  resulting  from  adverse  compressibility  effects  on  the 
advancing  blade. 

1.  Reducing  torsional  stiffness  and/or  torsional  inertia  decreases 
the  amplitude  of  the  vibratory  torsional  moments  associated  with 
retreating  blade  stall. 

2.  Moving  the  airfoil  elastic  and  center-of-gravity  axes  forward 

of  the  aerodynamic  center  axis  reduces  the  amplitude  of  the  high 
frequency  torsional  moments. 

3.  Increasing  blade  twist  generally  delays  the  airspeed  for  control 
load  buildup  but  increases  the  rate  of  control  load  buildup 
with  airspeed. 

i*.  Torsional  mode  shape  and  flatwise  and  edgewise  stiffnesses  do  not 
strongly  affect  high  frequency  torsional  moments. 

5.  An  increase  in  airfoil  static  stall  angle  reduces  vibratory 
pushrod  load.  Leading  and  trailing  edge  stall  airfoils  with 
the  same  stall  angle  result  in  approximately  the  same  level 
of  stall  flutter  response. 

With  regard  to  the  position  of  the  control  load  boundary,  the  results  of 
the  parametric  study  are  presented  graphically  in  Figure  87  where  the 
various  control  load  boundaries  are  compared  and  the  increase  in  allow¬ 
able  advance  ratio  with  each  independent  design  change  is  extracted  and 
plotted  versus  C^/a.  The  following  general  trends  are  evident: 

1.  Reducing  torsional  frequency  to  U . OP  by  reducing  blade  stiffness 
has  its  most  beneficial  effects  on  the  control  load  boundary  at 
high  C T/a. 

2.  Decresaing  torsional  inertia  while  maintaining  the  standard 
frequency  has  beneficial  effects  at  all  thrust  coefficients. 

The  maximum  increase  in  allowable  airspeed  of  i+0  knots  occurs 
at  the  intermediate  thrust  coefficient. 

3.  Forward  movement  of  the  blade  elastic  axis  and  increasing  blade 
twist  have  smaller  effects  on  the  pushrod  load  boundary.  High 
twist  decreases  the  stall  speed  at  low  C^/a  and  high  advance 
ratio. 
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REEVALUATION  OF  PRELIMINARY  STALL  BOUNDARY  PREDICTION  METHODS 


In  light  of  the  results  of  the  stall  boundary  generation  and  blade  design 
studies,  the  preliminary  stall  boundary  prediction  techniques  described 
earlier  were  reevaluated.  Tne  objectives  were  to  determine  (l)  whether 
any  of  the  simple  methods  can  be  expected  to  give  accurate  prediction  of 
the  stall  boundaries  of  an  aribtrary  rotor  design  and  (2)  whether  the  re¬ 
sults  of  the  parametric  design  study  can  be  used  to  extend  simple  design 
methods  by  including  the  effects  of  blade  structural  parameters  on  stall 
boundaries . 

It  is  clear  from  the  parametric  design  calculations  that  blade  structural 
design  changes  can  have  potentially  large  effects  on  stall  boundaries. 
Therefore,  boundary  prediction  techniques  which  rely  on  the  correspondence 
of  stall-related  flight  limitation  and  a  purely  aerodynamic  stall  criterion 
should  not  necessarily  be  expected  to  predict  accurately  the  difference  in 
operating  limits  for  two  rotors  with  widely  different  structural  design 
characteristics.  The  simple  aerodynamic  criteria  are  generally  based  on 
experience  with  a  given  type  of  rotor  and  become  less  valuable  as  departures 
from  past  rotor  designs  are  considered.  Also,  the  analytical  results  show 
that  many  types  of  stall  boundaries  exist  (control  loads,  blade  stress, 
power,  etc.)  and  that  design  changes  can  have  conflicting  effects  on  dif¬ 
ferent  boundaries.  Increased  blade  twist,  for  example,  expanded  the  CH-53A 
control  load  boundary  but  reduced  limiting  speeds  based  on  flatwise  stress. 
Information  such  as  this  is  not  likely  to  be  available  from  semiempirical 
design  approaches  which  are  generally  based  on  experience  with  one  type  of 
stall  boundary.  It  has  also  been  noted  that  prediction  of  the  effects  of 
blade  parameters  on  incipient  stall  boundaries  may  not  be  sufficient  for 
design  purposes.  For  example,  increasing  blade  twist  delayed  the  buildup 
of  pushrod  loads,  but  reduced  the  maximum  allowable  airspeed  at  low  values 
of  Cp/o  (Figure  7l). 

Because  it  is  desirable  to  have  rapid  but  accurate  means  of  estimating 
stall  boundaries  in  preliminary  design  trade-off  studies,  attempts  must  be 
made  to  integrate  the  results  of  the  flexible  blade  analyses  into  prelimi¬ 
nary  design  procedures.  The  effects  of  four  blade  parameters  on  the  stall 
boundaries  of  a  fairly  typical  low  twist,  symmetrical  airfoil  section  rotor 
have  been  evaluated  above.  Hopefully,  similar  analyses  can  be  performed 
which  extend  these  results  to  include  the  effects  of  other  parameters  such 
as  airfoil  camber,  tip  geometry,  and  control  system  characteristics.  If  a 
catalogue  of  such  results  was  available  to  the  blade  designer  the  effects 
of  rotor  structural  design  parameters  on  stall  boundaries  could  be  included 
in  preliminary  design  calculations.  One  of  the  various  aerodynamic  stall 
criteria  could  be  used  to  establish  a  basexine  stall  boundai  r,  and  then  cor¬ 
rections  to  the  boundary  positions  could  be  made  for  independent  structural 
design  changes.  For  example,  Ap  versus  Cp/cr  relationships  of  the  form 
shown  in  Figure  8j  could  be  used  to  modify  basic  stall  boundaries  to  account 
for  independent  blade  property  changes.  As  an  illustration.  Figure  88  shows 
the  measured  pushrod  load  limit  and  two  simply  calculated  stall  boundaries. 
The  first  curve  is  based  on  the  rotor  profile  torque  criterion  bCQ^/o  = 
0.008  and  the  second  on  ^  270°)  =  ^  methods  correspond  to 
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an  upper  stall  boundary  and  appear  to  give  good  agreement  with  the  vibra¬ 
tory  pushrod  load  structural  limit.  Also  shown  are  the  Ap  versus  Op /a  cor¬ 
rection  curves  for  fcur  parameter  changes.  Use  of  the  latest  analytic 
methods  to  generate  correction  curves  of  this  type  will  provide  significant¬ 
ly  increased  capability  for  estimating  stall  boundaries  by  allowing  the 
treatment  of  blade  parameters  not  covered  by  rigid  blade  analyses.  The 
complete  aeroelastic  analysis  can  then  be  used  to  verify  results  quantita¬ 
tively  and  to  investigate  the  effects  of  combining  two  or  more  design 
features . 
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CONCLUSIONS 


1.  Two-dimensional  stall  flutter  tests  indicate  that  reducing  blade  tor¬ 
sional  stiffness,  reducing  blade  torsional  inertia  and  moving  blade 
pitch  axis  forward  decrease  stall  flutter  induced  moments.  Inception 
of  stall  flutter  was  delayed  with  the  SC  1095  airfoil  relative  to  the 
NACA  0012  airfoil;  however,  once  initiated,  stall  flutter  loads  for 
the  two  airfoils  were  generally  comparable. 

2.  The  trends  shown  by  two-dimensional  wind  tunnel  stall  flutt  :•  tests 
can  be  predicted  using  two  unsteady  aerodynamic  models.  The  Time 
Delay  unsteady  aerodynamic  model  gives  a  conservative  prediction  of 
response  amplitude . 

3.  Acceptable  correlation  of  CH-53A  blade  stresses,  control  loads  and 
required  power  was  obtained  with  a  rotor  aeroelastic  analysis  employing 
variable  rotor  inflow  and  unsteady  aerodynamics.  Best  correlation  was 
achieved  using  the  a,  A,  B  unsteady  aerodynamic  model  of  Reference  1 

in  which  airfoil  unsteady  coefficients  are  expressed  as  functions 
of  angle  of  attack  and  its  first  two  time  derivatives. 

4.  Analytical  studies  to  determine  the  influence  of  design  parameters  on 
control  loads  generated  by  rotors  having  symmetrical  airfoil  sections 
generally  substantiated  two-dimensional  airi'oii  test  results.  That 
is,  reduced  torsional  blade  stiffness,  reduced  torsional  inertia, 
increased  airfoil  static  stall  angle  and  forward  movement  of  the  blade 
pitch  axis  all  led  to  reduced  dynamic  control  loads. 

5.  Increasing  blade  twist  delayed  the  control  load  stall  speed  at  high 
Gp/o  but  increased  the  rate  of  control  load  buildup  with  airspeed. 
Increasing  twist  increased  blade  flatwise  moments  substantially. 

6.  Current  techniques  for  predicting  stall  boundaries  which  are  based  on 
aerodynamic  stall  parameters  alone  are  not  equally  applicable  to 
blades  having  differing  dynamic  properties. 
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RECOMMENDATIONS 


1.  Further  analytical  studies  should  be  conducted  to  examine  the  combined 
effects  of  the  design  parameters  which  have  been  found  to  be  promising 
in  this  study.  In  addition,  other  design  parameters  not  considered  In 
this  study,  such  as  airfoil  camber,  blade  tip  sweep  and  number  of 
blades,  should  be  evaluated. 

2.  Model  rotor  tests  should  be  conducted  to  substantiate  the  predicted 
beneficial  effects  of  torsional  stiffness,  torsional  inertia,  twist 
and  airfoil  contour  on  stall  boundaries. 

3.  Further  work  should  be  done  to  refine  the  Time  Delay  unsteady  aerody¬ 
namic  model  or  on  similar  techniques  which  do  not  require  test  data 
for  the  airfoil  being  modeled. 

4.  The  sources  of  blade  torsional  moments  experienced  at  V  =  180  deg 
should  be  studied.  In  this  area,  high  Mach  number  unsteady  airfoil 
pitching  moment  data  should  be  obtained  and  improved  procedures  tor 
predicting  the  lifts  and  moments  associated  with  close  blade-cortex 
passage  should  be  developed.  The  influence  of  spanwise  velocity 
components  should  al;o  be  assessed. 

5.  The  present  aeroelasiio  analysis  should  be  extended  to  include  the 
effects  of  fixed  control  system  dynamics  on  blade  response.  Tests 
involving  changes  to  fixed  system  design  should  be  performed. 
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TABLE  I.  TWO-DIMENSIONAL  AIRFOIL  TEST  POINTS 


IrKF  *  0.0005  slug  ft  /ft t  cT  =  8  deg 


Airfoil 


HACA  0012 
and 

SC  1095 


Torsional 

Natural 

Frequency 

“’e  •  CFS 


Oscillatory 
Frequency 
ft  ,  cps 


Torsional 

Inertia 

Ratio 


Mean  Incidence 
Angle 

V  de« 


7* 

11* 

lit* 

18 

7 

11 

lU 

18 

7 

11 

lit 

18 

7* 

11* 

lit* 

18 

7 

11 

lU 

18 

7 

11 

1U 

18 

7* 

11* 

lU* 

18 

7 

11 

lit 

18 

7 

11 

lU 

18 

7* 

11* 

lU* 

18 

7 

11 

lit 

18 

7 

11 

lit 

18 

7* 

11» 

lit* 

16 

7 

11 

lU 

18 

7 

11 

lU 

18 

7* 

11* 

lit* 

18 

1  7 

11 

lit 

18 

!  7 

11 

lit 

18 

j  T« 

11* 

lU* 

18 

7 

11 

lit 

18 

7 

11 

lit 

18 

7* 

11* 

lit* 

18 

•Analytic  Cases 


TABLE  II.  MODEL  AIRFOIL  ELASTIC  RESPONSE  AMPLITUDES 


NACA  0012  Airfoil  Si  =  10  cps 

.  xpA  =  0.25 

Torsional 
Inertia  Ratio 

I/I:REF 

Torsional 
Natural  Freq. 
W0.  CP"» 

Mean  Incidence 
Angle  om.  deg 

A6. 

deg 

6 

HPTP 

deg 

L.O 

50.7 

L.O 

50.7 

L.O 

50.7 

L.O 

50.7 

L.O 

70.0 

L.O 

70.0 

L.O 

70.0 

L.O 

70.0 

L.5 

5^.8 

L.5 

5^.8 

L.5 

51+.8 

L.5 

51*.8 

L.5 

69.0 

L.5 

69.0 

L.5 

69.0 

L.5 

69.0 

NACA  0012  Airfoil  Si  =  10  cps ,  xpA  =  0.22 


«  Stall  Flutter  not  present 


TABLE  HI.  MODEL  AIRFOIL  ELASTIC  RESPONSE  AMPLITUDES 

CALCULATED  USING  a.A.B  UNSTEADY  AERODYNAMICS 


NACA  0012  Airfoil  0  =  10  cps, 

,  xpA  =  0.25 

Torsional 

Torsional 

Mean  Incidence 

A6  , 

e 

Inertia  Ratio 

Natural  Freq. 

Angle  a^,  deg 

deg 

JjPTP 

I/:IREF 

u>6,  cpra 

deg 

NACA  0012  Airfoil  Q  *  10  ops,  xpA  =0.22 


0.32 

0.87 

0.67 

1.20 

l.lU 

1.73 

O.lU 

0.30 

o.uu 

0.85 

0.58 

1.00 

0.2U 

0.65 

0.50 

0.93 

0.71* 

1.15 

0.10 

0.22 

0.3U 

0.59 

0.35 

0.60 

TABLE  TV.  MODEL  AIRFOIL  ELASTIC  RESPONSE  AMPLITUDES 
CALCULATED  USING  TIME  DELAY  AERODYNAMICS 


NACA  0012  Airfoil  ft  =  10  cps,  xpft  =  0.25 


Torsional 
Inertia  Ratio 

1 7  XREF 

Torsional 
Natural  Freq. 
gj0,  cpm 

Mean  Incidence 
Angle  <*M,  deg 

deg 

0 

LPTP 

deg 

1.0 

50.0 

7 

0.  It  It 

0.82 

1.0 

50.0 

11 

1.28 

2.11 

1.0 

50.0 

lit 

2.05 

l.olt 

1.0 

70.0 

7 

0.22 

O.ltP 

1.0 

70.0 

11 

0.86 

l.ltO 

1.0 

70.0 

lit 

1.02 

20 

r**' 

<— 1 

1.5 

50.0 

7 

0.19 

0.7? 

1.5 

50.0 

il 

1.0 

1 . 36 

1.5 

50.0 

lit 

1.68 

1 

1.5 

70.0 

7 

0.18 

0.36 

1.5 

70.0 

11 

0.65 

0 . 80 

1.5 

70.0 

lh 

0.70 

0.9? 

NACA  0012  Airfoil,  ft  =  10  cps,  xt 


1.0 

50.0 

7 

0.32 

0 .  .5? 

1.0 

50.0 

11 

0.85 

1.57 

1.0 

50.0 

lit 

1 .66 

;>.'.8 

1.0 

70.0 

7 

LT\ 

< — 1 

O 

0.15 

1.0 

70.0 

ll 

0.7? 

1  .'5i 

1.0 

70.0 

lit 

0.05 

1.23 

1.5 

50.0 

7 

0 . 2  0 

0 . '* 

1.5 

50.0 

u 

1 . 1  !i 

1.20 

1.5 

50.0 

lit 

1 .60 

2.13 

1.5 

70.0 

7 

0.15 

0.15 

1.5 

70.0 

11 

n.57 

.7  it 

1.5 

70.0 

1^ 

O.f  5 

0.87 

TABLE  IV  -  Concluded 


X  1095  Airfoil  52  =  10  cps ,  x  =  0.05 


Torsional 


Torsional 
Natural  Freq. 
wQ ,  cpm 

Mean  incidence 
Angle  aM,  deg 

t  o 

:  v 

deg 

o 

H.PTP 

deg 

50.0 

7 

o.i>3 

0.1,3 

50.0 

11 

Q  .  )(  J 

0 . 1*  3 

50.0 

Id 

O 

Ln 

3.0  It 

70.0 

7 

0.1  ' 

0.19 

70 .0 

11 

0 . 1  1 

0.19 

70.0 

Ik 

1 .07 

1 . 81. 

50.0 

7 

'  \  j’  7 
•  ■  1 

p .  37 

50.0 

11 

0.  -7 

0 . 37 

50.0 

la 

1 .  J8 

2.1.8 

70 . 0 

7 

0.17 

0.17 

70.0 

11 

0.  i? 

0.17 

70.0 

Ik 

0.73 

1  .]■. 

SC  1095  Airfoil  .9  =  1 0  crs,  x  =  ).9i 


ALL  DIMENSIONS  IN  INCHES 


TOP  VIEW  WITH  NACA  0012  FULL  SPAN  AIRFOIL  INSTALLED 


SIDE  VIEW 


Figure  1.  Acoustics  Research  Wind  Tunnel  Test  Section. 
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Figure  11.  Airfoil  Steady-State  Drag  Characteristic 
Corrected  for  Open  Jet  Tunnel  Effects. 


;;ACA  J012  Model  Airfoil  Angle  of  Attack 
Tin<=-  Histories. 


Y.xi' 


=  Iref 


1  gur< 


T-  counting  begins  t 


=  moment  stall  time  constant  r  =  2.0 

.  s 

.  C.P.  begins  to  move  rearward  with  time  •  ward  C.P 
.  c  „„  is  eliminated 
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Increased  by  a  factor  of  3 
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T‘°nel  =  T?n  *  2AV“°n) 
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Figure  l3.  rime  DeLay  Unsteady  Aerodynamic 
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NACA  0012  AIRFOIL 


IV ’ 1  Iffff 


orrelation  of  JACA  0012  Airf  :i  3  spouse  Ar.flitude  Obtained 
'ith  Two  Unsteady  Aerodynamic  i*t'n:ds. 
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Selected  Measured  Stall  Boundaries 
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he  Buildup  of  CH  -53A  Pushrod  Loads  With  Airspeed. 


VIBRATORY  EDGEWISE  VIBRATORY  FLATWISE 

STRESS,  ±  PSI  STRESS,  ±  PSI 


0  .2  A  .6  .8  1.0 


NONDIMENSION AL  RADIAL  STATION 


Figure  UO.  Effect  of  Nonunifonn  Inflow  on  Correlation  of  the  Radial 

Distribution  of  CH-53A  Blade  Stress  at  137  Kn  and  1*2, 00C  Lb 
Gross  Weight. 
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00  Lb  Gross  Weight. 
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Blade 


relation  of  CH-53A  Blade  Stress  With  a.  A,  B  and  Time  Delay 
teady  Aerodynamics  at  137  Kn  and  U2»000  Lb  Gross  Weight. 


I  ROD  LOAD, 


MAXIMUM  VIBRATORY  EDGEWISE  MAXIMUM  VIBRATORY  FLATWISE 
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Figure  U8.  Correlation  of  CH-53A  Vibratory  Blade  Stress  at 
42,000  Lb  Gross  Weight. 
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-the  Buildup  of  CH-5^B  Pushrod.  Loads  With  Airspeed. 


2000 
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(N0ISN31  +  )  81  QV01  OOMHSnd 
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Figure  51.  Correlation  of  CH-S^B  Pushrod  Load  Tine  Histories  at  115  Kn 
and  U7.000  Lb  Gross  Weight. 


STALL  FLUTTER  SPIKE  AMPLITUDE.  VIBRATORY  PUSHROD  LOAD. 


Figure  52.  Correlation  of  Measured  and  Calculated  CH-5^B  Pushrod 
Load  Amplitudes  at  U7.000  Lb  Gross  Weight. 
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MAXIMUM  VIBRATOR  V  MAXIMUM  VIBRATORY 

EDGEWISE  STRESS  ,±PSi  CLATWISE  STRESS,  ±PS 


MAXIMUM 


Figure  ,  Flight  Conditions  Used  To  Generate 
CH-53A  Ctall  Boundaries. 
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Figure  55.  The  Buildup  of  Calculated  Pushrod  Loads  With  Airspeed 
for  the  CH-53A  at  Gp/o  =  0.083. 
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PUSHROD  LOAD,  LB  (+TENSION) 


-1000 


Figure  56.  The  Buildup  of  Calculated  Pushrod  Loads  With  Thrust 
Coefficient  at  u  =  0.33. 
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VIBRATORY  pushrod  load,  ±  lb 


Figure  57-  Agreement  Between  Test  and  Analytical  Pushrod  Load  Amplitudes 
and  Dtall  Boundaries  for  the  CH-53A. 
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■'  igure  59-  Calculated  Buildup  of  CH-53A  Vitratory  Blade  Stres 
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Figure  6l.  Calculated  Stall  Boundaries  for  the  CH-53A  Rotor. 
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Figure  67.  Effect  of  Forward  Movement  of  Elastic  and  Center- 
of-Gravity  Axes  on  CH-53A  Vibratory  Pushrod  Loads 
and  Required  Power. 
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Figure  i/>.  Effect  of  Forward  Movement  of  Elastic  ana  Center- 
o  f-'Jravi  ty  Axes  on  Cli—  S  3A  Stall  Boundaries. 
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2.000  LB 


Effect  of  Increased  Blade  Twist  on  CH-53A  Pushrod  Load  Time  Histories 


MAIN  ROTOR  HORSEPOWER  VIBRATORY  PUSHROD  LOAD,  *  LB 


Figure  71.  Effect  of  Increased  Blade  Twist  on  CH-53A  Vibratory 
Pushrod  Load  and  Required  Power. 
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Figure  72.  Effect  of  Increased  Blade  Twist  on  CH-53A  Blade 
Stresses  and  Available  Control  Moaent . 
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Figure  73.  Effect  of  Increased  Blade  Twist  on  CH-53A  Stall  Boundaries. 
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Figure  7 6.  Variation  of  Torsional  Response  Frequency 
With  Natural  Frequency. 
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ROTOR  HORSEPOWER 


Figure  79.  Effect  of  Reducing  Torsional  Natural  Frequency  on 
CH-53A  Blade  Stresses  and  Available  Control  Moment. 
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STALL  FLUTTER  SPIKE  AMPLITUDE, 


Figure  83.  Trends  of  Stall  Flutter  Spike  Amplitude  With  Torsional 
Frequency  and  Stiffhess  Ratios. 
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Figure  85.  Effect  of  Decreased  Torsional  Inertia  and  Torsional 
Stiffness  on  CH-53A  Blade  Stress  and  Available 
Control  Moment. 
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VAXIMUM  THRUST  COEFFICIENT  —  SOLIDITY  RATIO .  CT  /a\  Mflx 
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Figure  87.  Effect  of  Blade  Design  Changes  on  the  CH-53A 
Control  Load  Boundary. 
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MAXIMUM  ThRUjT  „DEFFiC  IENT 
SOLIDITT  RATIO  ,Cr/cr\ 

'  I  MAX 


CHANGE  IN  CONTROL  LOAD  LIMIT  ADVANCE  RATIO,  Ag. 


Figure  85.  Modifications  to  the  Control  Load  Stall  Boundary  To  Be  Used  in 
Preliminary  Design. 


LIST  OF  SYMBOLS 


A  dimensionless  angular  velocity 

b  number  of  blades 


B  dimensionless  angular  acceleration 

B^g  longitudinal  cyclic  pitch,  deg,  positive  for  increased  pitch 

angle  at  270  deg  azimuth 

c  airfoil  chord,  ft 

C  mechanical  damping  per  unit  span,  (ft-lb-sec/rad)/ft , 

Equation  (5) 

c^  section  drag  coefficient 

c^  section  lift  coefficient 

clP0T  potential  flow  section  lift  coefficient.  Figure  l8 

C1SEP  fully  separated  section  lift  coefficient.  Figure  18 

cm  section  pitching  moment  coefficient 

cmc/l+  section  pitching  moment  about  quarter  chord  position 

cmxpA  section  pitching  moment  coefficient  around  pitch  axis 

C.P.  airfoil  section  center  of  pressure.  Figure  18 

^Q,D  rotor  profile  drag  torque  coefficient 

T 

CT  rotor  thrust  coefficient  7Tp;2p  (ar)2 

ci,  C2  constants.  Equation  13 
f  frequency,  cycles/sec 

GJ  blade  torsional  stiffness  per  unit  length,  in. -lb/rad 

h  Jet  height,  ft.  Equation  (l) 

hD  density  altitude,  ft 

I  airfoil  torsional  moment  of  inertia  in  pitch  per  unit  span, 

(ft-lb-sec2)/ft,  Equation  (5) 
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reference  value  of  airfoil  torsional  moment  of  inertia  per  unit 
span, (ft-lb-sec2)/ft 

airfoil  or  blade  torsional  inertia,  ft-lb-sec2 
reduced  frequency 

pitch  spring  stiffness  per  unit  span,  (ft-lb/rad)/ft , 

Equation  (5) 

equivalent  spring  stiffness  of  flexure,  in. -lb/rad 
load  factor.  Equation  (10) 

Mach  number 

applied  aerodynamic  moment  per  unit  span,  ft-lb/ft,  Equation  (5) 

nondimensionalized  blade  aerodynamic  pitching  moment,  positive 
for  noseup  twisting  moment,  Equation  ( 11 ) 

nondimensionalized  blade  center-of-gravity  offset  torsional 
moment,  positive  for  noseup  twisting  moment,  Equation  (ll) 

nondimensionalized  blade  torsional  inertia  moment,  positive  for 
noseup  twisting  moment,  Equation  (ll) 

nondimensionalized  blade  structural  moment,  positive  for  noseup 
twisting  moment.  Equation  (ll) 

nondimensionalized  blade  torsional  moment  caused  by  shear  de¬ 
flection  effects,  positive  for  noseup  twisting  moment,  Equation 
(11) 

flexure  torsional  moment,  in. -lb 

normal  rotor  rotational  speed,  rpm 

period  of  model  airfoil  driving  motion,  sec 

rotor  pitching  moment,  ft-lb,  positive  for  nosedown  moment 

rotor  section  radial  coordinate,  ft 

rotor  radius ,  ft 

time,  sec 

rotor  thrust,  lb 

applied  torque  per  unit  span,  ft-lb/ft,  Equation  (5) 
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Uq  free-stream  velocity,  ft /sec 

Vc  calibrated  airspeed,  kn 

W  aircraft  gross  weight,  lb.  Equation  (10) 

^AC  nondimensional  chordwise  aerodynamic  center  location 

Fcq  nondimensional  chordwise  center-of-gravity  location 

**EA  nondimensional  chordwise  elastic  axis  location 

xpa  nondimensional  chordwise  pitch  axis  location 


aM  airfoil  mean  angle  of  attack,  deg  or  rad 

cT  angular  amplitude ,  deg  or  rad 

static  pitching  moment  stall  angle,  deg  or  rad.  Figure  18 

o^,  angle  of  attack  above  which  center  of  pressure  is  assumed 

fixed,  deg  or  rad.  Figure  l8 

5ci  C1  correct,ion  f°r  scaling  unsteady  lift  coefficient.  Equation  (7) 

6c  c  correction  for  scaling  unsteady  pitching  moment  coefficient, 

m  Equation  (8) 

6m  pitching  moment  coefficient  correction  curve  for  leading  edge 

stall  airfoil, rad.  Equation  (12) 

6a^  angle  of  attack  correction  for  scaling  unsteady  lift  coefficient, 

rad,  Eauation  (j) 

6c.  angle  of  attack  correction  for  scaling  unsteady  moment  coefficient, 

rad,  Equation  (8) 

Aa  angle  of  attack  correction  for  open  jet  effects,  rad.  Equation  (l) 

Ac^  open  Jet  drag  coefficient  correction.  Equation  (2) 

Acmc/l|  open  Jet  pitching  moment  correction,  Equation  (3) 

AtgEp  time  between  static  stall  and  dynamic  lift  stall,  sec.  Figure  18 

AO^  half  amplitude  of  airfoil  torsional  deflection,  deg  or  rad. 

Figure  ih 

0  airfoil  or  rotor  torsional  elastic  deformation,  deg.  Equation  (k) 
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blade  linear  twist  angle,  deg,  positive  for  tip  twisted  up 
relative  to  root 

half  amplitude  of  overall  airfoil  torsional  deflection,  deg 
cr  rad,  Figure  1^ 

rotor  advance  ratio 

two-dimensional  aerodynamic  damping  in  pitch 

O 

air  density,  slugs/ft 

rotor  solidity 

nondimensional  time,  t/P 

nondimens ional  elapsed  time.  Figure  l8 

nondimensional  time  delay  constant  for  lift  stall 

nondimensional  time  delay  constant  for  moment  stall 

blade  root  azimuth  angle  measured  from  downstream  blade  position 

torsional  natural  frequency,  cycles/sec 

torsional  natural  frequency  ratio,  u>e/ft 

airfoil  or  rotor  driving  frequency,  cycles/sec 

first  derivative  with  respect  to  time 

second  derivative  with  respect  to  time 

first  generalized  torsional  mode 


